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FOREWORD 


This document constitutes a final report on the application of the 
blast-wave theory of meteoroid impact to the problem of space radiator 
design. Other publications entitled "Comments on the Solution of the Spall- 
Fracture Problem in the Approximation of Linear Elasticity", CAL Report 
AI-1821-A-3, NASA CR-54250, and "On the Possibility of Simulating 
Meteoroid Impact by the Use of Lasers", CAL Report AI-1821-A-1, NASA 
CR-54029, have previously been generated under this NASA-sponsored pro- 
g ram. 

The present report is mainly concerned with the phenomena of shock 
propagation and fluid flow due to impact. The application of these results 
to the prediction of the deformation remaining after the fluid motion ceases 
is taken up in the report on spall fracture, mentioned above. 
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ABSTRACT 


This report presents analytical solutions for the£propagation of 
shock waves in dense targets struck by hypervelocity projectiles. The 
solutions considered deal mainly with the high-pressure phases of defor- 
mation, during which the material strength can be neglected. Several 
methods of solution are described. Their advantages and limitations are 
discussed, and their accuracy is estimated by comparing their predic- 
tions with experiment and with the results of numerical solutions. To 
illustrate the application of these approximate theories, an example is 
given of shock propagation in porous targets. 
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INTRODUCTION 


The deformation of a solid target by the impact of a high-speed 
projectile is an extremely complex phenomenon. Attempts to predict the 
behavior of real materials during such an intense disturbance are com- 
plicated by the presence of a variety of nonelastic, nonlinear effects. 

In recent years, however, interest in impact at the extreme velocities 
of meteoric particles (10-70 km/sec) has made it possible to neglect 
material strength during much of the deformation. This approximation 
is suggested by the fact that the pressures generated at impact are typi- 
cally measured in megabars, and far exceed the strength of any common 
material. While the applicable equations of inviscid motion are consid- 
erably simpler than those in which material strength is accounted for, 
nevertheless their solution presents a formidable problem in time- 

dependent, multidimensional, compressible fluid flow. Numerical solu- 

1 2 

tions of such problems for specific cases have been reported by Bjork 

3 4 

and by Walsh and his coworkers. ’ The latter results, in particular, 
have revealed the existence of certain scaling laws, but extensive corre- 
lations of the numerical results have not yet appeared. Parallel analytical 
treatments, while less specific, have also identified these scaling laws, 

5 

and have been used to correlate results from a variety of cases. The 
main purpose of this report is to review several of the more promising 
analytical approaches that have been used, to present some extensions of 
them, and to show their range of application. 

The solutions described below deal exclusively with the inviscid- 
fluid model; the strength mechanisms by which the target material is 
eventually brought to rest are omitted. Consequently, the solutions offer 
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no direct information about target damage, but describe only the shock 
propagation and the flow behind the shock. Some implications of these 
results on the question of target damage are presented in the Concluding 
Rema rks. 

One of the main obstacles to finding analytic solutions for shock 

propagation in solids is that the powerful method of similarity cannot be 
5 

used. The reasons for this fact, which is associated with the equation 
of state of the material, are briefly reviewed in Section I. Various forms 
of the state equation are considered in Section II, and the particular form 
used in the present report is described in detail. The three types of non- 
similar solutions considered are taken up in Section III. The first of 
these, developed by Oshima 5 for a perfect gas, has perhaps the greatest 
analytical justification, and is here extended to the case of solids. 

The second is the "Quasi-Steady" solution used previously 5 to 
achieve a fairly accurate correlation of impact-generated shock waves. 
This theory is also extended, so as to apply to the weaker, later stages 
of shock propagation, and is compared extensively with experiment and 

with numerical solutions. The third method reviewed was developed ori- 
7 

ginally by Porzel for the case of atmospheric blasts, and was later 

g 

extended by Zaker to the problem of a nuclear explosion inside a solid. 
The Po rzel-Zaker theory is rederived here, in the light of new informa- 
tion now available about the equation of state of solids. 

The application of these theories is illustrated in Section IV by 
considering shock propagation in porous targets, a situation which appar- 
ently has not yet been treated numerically. The report concludes with 
rema rks intended to indicate what can and cannot be done at present by 
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analytical means, and to call attention to several lines of research that 
may be fruitful for achieving an improved analytical understanding of 
deformation due to hypervelocity impact. 
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I. SIMILARITY CONSIDERATIONS 


When material strength is neglected, the equations expressing the 
conservation of mass, momentum, and energy reduce to the Euler equa- 
tions of an inviscid, compressible fluid. Assuming that all motion is 
symmetric about an axis normal to the target surface, and using a spherical 
coordinate system with its center at the impact point, the equations of 
motion are 



+ U& + 

dr r 9B 



J_ &£ _i_ T>U_ 

r s>e f 


■h ~~~ — O 


( 1 ) 


da 

at 



ur £g 

r ae 



/ ar 


— o 


( 2 ) 



Here U. and uf denote velocity components in the f - and Q -direc- 
tions, ^ denotes the pressure, p the density, and e the internal 
energy per unit mass. The index 7 ) appearing in the continuity equation 
must, strictly speaking, be set equal to 2 , since these equations apply 
only for a spherical problem. Its value is left arbitrary for the moment, 
since, when the Q -variation is discarded (as will be done shortly) the 
correct forms for the cylindrical and planar cases can be recovered by 
setting - 2 ) = 1 or 0 . 
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These equations are completed by the equation of state, which is 
taken for the moment in the form 

e = < 5 > 

At the shock, the Rankine -Hugoniot conditions apply. Thus 


U, 


f',A 




fo, hj e o 

U-o — O 




( 6 ) 


'f' t° ~ A ^3 cL , 


(7) 


- ~ 



( 8 ) 


where ££ is the velocity of any portion of the shock normal to itself, and 
subscripts o and 1 denote conditions ahead of and behind the shock. 

The analysis presented in this report is based on the assumption that 
the shock is always a hemisphere, with center at the impact point on the 
surface. Such an assumption is supported by photographs of impact- 
generated shock waves in transparent targets, ^ by internal meas- 

12 13 

urements of shock shape in wax targets, * and by computer solutions 
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for impact into solids, at a sufficiently long time after impact. ^ 

For impact into a gaseous target, the results of Walsh et al 4 indicate a 
spherical shock whose center lies inside the original target surface, at 
about one-third the distance traveled by the shock along the axis of sym- 
metry. 

The problem of predicting the flow field behind the shock wave as 
it propagates into the target requires the solution of Eqs. (1) - (4), which 
contain two spatial coordinates and time. One of the most fruitful 
approaches to such a complicated problem has been the method of simi- 
larity, which assumes that properly normalized distributions of the 
various physical quantities (pressure, density, etc. ) at each instant are 
the same when viewed on a scale defined by the shock radius at that in- 
stant. Thus each quantity, instead of varying independently with time and 
with the spatial coordinate C (distance from the origin), is assumed to 
be a function only of the similarity variable. 

1 = */&,(*) (9) 

where R s (i) is the instantaneous distance of the shock from the origin. 
Mathematically speaking, the objective of the similarity assumption is to 
suppress time as an explicit independent variable. To see how this is 
done, and to set the stage for the introduction of Oshima's method below, 
it is instructive to change independent coordinates from f , ■& £ 

to y , & , g s (£) , and to redefine the dependent variables in terms 

of certain dimensionless functions: 
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( 10 ) 


U- = 4 Hi, f =/£ 4V(Y e 

^ <i) /> =/. '/•’('% e, ^) 

£-e„ = 4 ^ 61 ,^^) 

In these terms, the basic differential equations become 




£j 

^ ae 


fi ^ + ?>-£ -t- co-& e) = 








( 11 ) 


4" 


j> + (V 


•n 




1 a© 


aj*~ . 1 ^ _ ip 5^ 

? + fS[--^3e, 


( 12 ) 



The assumption of similarity amounts to neglecting the right-hand sides 
of these equations -- i. e. , it is assumed that the functions , Od > 
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£ » • an< l ^ do not change with the scale of the disturbance. 

• • 

In addition, similarity requires that the combination be a 

constant, which is satisfied if the shock radius is proportional to a power 
of the time 

rj 




(15) 


The only remaining source of time dependence then comes from the 
boundary conditions. The Rankine -Hugoniot conditions, written in terms 
of the dimensionless functions introduced above, become 


<4 (>-&) 



where conditions ahead of and behind the shock are denoted by subscripts 

O and 1 , respectively. These equations show that an explicit time 

dependence can enter in either of two ways. The first occurs whenever 

the pressure ahead of the shock is comparable to , which is of 

the order of the pressure behind the shock. The second possibility of 

time dependence comes from the state equation. As pointed out, for ex- 
14 15 

ample by Kynch and by Sedov, similarity solutions are possible only 
if the equation of state has the form 


e =• 'f <f (f>) 


AI-1821 -A-2 


8 



where ^tj^) an y ^ unct i° n of the density. When the state equation has 

this form, and when the term is negligible, it is possible to 

eliminate 


% 


between the last two of Eqs. (16), and the result is an 


expression for the density ratio across the shock , which remains 


constant 


VCf'M) - 1 + ir 

The equation of state of a perfect gas 

e = -t 


-= o 


08) 


(y-'Y 


(19) 


is of the form of Eq. (17). Thus, departures from similarity in the 
perfect-gas problem arise only from the term * The situa- 

tion is just the opposite for shock propagation in solids. There, the term 
jj? 1 * is certainly negligible whenever the fluid-mechanical model 

itself is appropriate. On the other hand, the state equation appropriate 
to a high-pressure solid does not match the form required byEq. (17), 
except in the limit of extremely high pressures, where, as will be seen 
below, the state equation approximates the required form, and the density 
ratio across the shock approaches a constant value. It is a matter of 
common experience, however, that the density ratios produced by impact- 
generated shock waves are far from constant. For impact speeds on the 
order of 1 0 km/sec, they may range as high as 2 or 3 at the impact point, 
but they drop rapidly as the wave expands, and quickly approach the 
elastic-wave limit of 1. Only the very early stages of impact at very high 
speed experience the constant density ratio required for similarity. Thus, 
as pointed out in Ref. 5, impact-generated shock propagation in solids is 
characteristically nonsimilar. 
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II. EQUATION OF STATE 


The nonsimilar behavior of shock propagation in solids is caused by 
the form of the state equation. For the range of pressures encountered 
in hypervelocity impact, the Mie -Griineisen equation of state is commonly 
used. In most numerical solutions, this equation is expressed in consid- 
erable detail, requiring the specification of a large number of constants. 
For present analytical purposes, a less cumbersome expression is desir- 
a ble, this section describes an approximate form of the Miie-Gruneisen 
equation used in the analyses that follow. A complete description is given 
of the shock-wave relations and thermodynamic properties for such a 
state equation. 

A. General Form 


The Mie -Griineisen equation is usually written in the form (see Ref. 


16, Eq. (33)) 

(/>) 


P r (f) 


( 20 ) 


where the subscript c denotes the cohesive contribution, and where 
r(fi is the Griineisen factor, which depends weaklyon . It is con 

venient to express the cohesive quantities in terms of measured shock- 
wave data; along the Hugoniot, Eq. (20) takes the form: 


/"/>) 


( 21 ) 


Subtracting the left and right sides of this from the corresponding sides 
of (20) gives 
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e - e H U) - -i — inJ^l 


(22) 


This can be rearranged in the form 



where 


(23) 



(24) 


The presence of the second term in Eq. (23) renders a similarity solution 
inappropriate (compare with Eq. (17)) except in the limit of extremely 
high pressure, where A is ne gli§ikle compared to the first term. 

B. Shock-Wave Relations for General f* , A 

The Rankine-Hugoniot condition expressing conservation of energy 
across a shock is 


P _ p — 

i c o — 



(25) 


If all energies are measured with respect to , the equation of state 

may be rewritten 

— L(f) (26) 

Equating these two expressions gives the Hugoniot curve in the form 


e-e, = 


L 


/ r 0) 



(27) 
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. An alternate 


where the approximation has been made that 


<r<r 


form can be derived by using the other two Rankine-Hugoniot conditions 






(28) 


This leads to 


{' 


zA 


)(/■/ Irtfi J/„ ( r( r .) 


I + 


= o 


(29) 


This relation is shown in Fig. 1. If information about the functions r(/>,) 
and A (f,) is available (perhaps from a theoretical estimate of the 
cohesive terms) this figure can be used, for a given value of />, , to 

find the corresponding value of . 

C. Specialization for a £__ L _£ Medium 

The steps indicated above can be reversed, and the measured 
Hugoniot data can be used to investigate the quantities r(^) and • 

In particular, if l ~ is known, Eq. (24) gives /\ (/)) in terms of the 
Hugoniot data. These relationships take on an especially simple form 
for a substance whose Hugoniot displays a linear relation between the 
shock speed and the particle speed behind the shock 


U s - <2 +■ SUL, 


(30) 


num- 


Such a substance is here referred to as a C , 5 medium. A large 
ber of materials are well approximated by such a relation. Reference 16, 
for example, presents a collection of values of <2 and S for various 
materials. 


As pointed out, for example, by Tillotson. 
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for such a material is 


In what follows, the explicit form of 
derived. From Eq, (24) 




(31) 


For the moment, is assumed to be a known constant; its precise value 

will be specified later. For a C , S medium 


lft4 ~ fo 4s 

Thus 



[l- s + s h/f> ] 


(32) 


2[s-(s-i)/^ 0 ] 1 ' 


(33) 


In the limit of extremely strong shocks, where — *-<39 , the re- 

sults of Fig. 1 indicate that the density ratio should approach + 1. 
On the other hand, for a C ,S medium, 


<09 

The constant value to be used for f is chosen so as to make these two 
limiting values for the density ratio equal, i. e. , 

■jr + -I - , or r- 2^s-0 , s = / + £> (34) 

It would be more accurate to allow r ( to vary from this value at 
the limiting density ratio to the value 25 - 1 at normal density, ^ per- 
haps in a linear fashion 


A Cs-.)“*4 + 


S— I 




c. 
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(35) 


r(f)= zs-i _ (s-i)(£-i) 

but for purposes of convenience, the constant value f~ = 2S -2 has been 
used throughout the analysis that follows. 

With this specification, the final form of the state equation used 
below is 


C* 



/ 

2 (s —\ ) 


r° 

s-ts-OS/fi 


(36) 


This representation has the advantage of generality. It shows that the 
equation-of-state data for a large number of materials can be correlated 
in terms of the dimensionless pressure U.o , density A/ A 0 , and 
specific internal energy — , with only S remaining as a pa ram- 

jjc 

eter. Figure 2 is a plot of this equation for the cases S = 2 and 
S = 1. 5. The qualitative resemblance to the specific state equation of iron 
used in Bjork's calculations (see Fig. 1 of Ref. 18) is obvious. Figure 2 
shows the Hugoniot, which can be found from 


e H -e 0 _ , 
c 2 - ^ 


(±1 

[s-rs-oA ^] 2 





(37) 


as well as a number of isentropes, about which more will be said below. 


The fact that the internal energy scales with C z can be used to esti- 
mate how much the Hugoniot is changed by heating a target, such as a 
heated radiator tube, for example. Changes in e o due to heating are 
given by C^AT . Even for a temperature rise of 1000°K, C V AT / c *- 
is on the order of a few percent. Thus the shock-wave properties of a 
given material are not appreciably altered by heating it, and it might be 
expected that changes in impact damage in heated targets are due mainly 
to the effects of temperature on strength. 
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The Hugoniot conditions can also be written conveniently in terms of the 
shock Mach number Ml — • A few values of Ml are shown on Fig. 

2. The corresponding values of pressure, density, internal energy, and 
particle velocity are 


A = — 


= (j±^L 

Z 5 l 


fn _ ) 

^ 5 


(40) 


u, __ M- / 

<=• £ 


(41) 


The representation afforded by Eq. (36) ceases to be valid in the 
limit of extremely high pressures, where the linear , U ^ relation 

no longer holds. As noted below, Eq. (36) is also incorrect at low pres- 
sures. It is, however, a simple and useful approximation in the inter- 
mediate range of pressures -- i. e. , in the range where the approximation 
of Eq. (30) is appropriate. 

D. Isentropes; Low-Pressure Behavior 
Along an isentrope, the relation 

^-=4- (42) 

f 
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together with the general form of the Mie -Gruneisen relation, gives the 
following linear differential equation 


— (e-e*) - 

whose solution is 


(e-e>) = 


= const +. r ^ A ^ ~ 


For the special case of a C , 3 medium, taking /~ = 2 ( S - 1) and 


using 


^--1 

AgO = _/ A 1 

C v Z(s- 1) s- 

the equation for the isentropes is found to be 


^ = U 
(a 


z( S-l) r ^ 

) /fi + ;s ) 


where 




Equation (46) defines the isentrope passing through the point — g ° — £ 

a c- 

~ constant / can be evaluated in terms of conditions along 

the Hugoniot; an isentrope starting from the point ^ will 

have 

/v , K-z(s-f) / 

4 / M e*zs- \ ,48, 
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The function iD can be evaluated explicitly whenever 25-1 is 


an integer. In particular 


D ^ il,r)= tty. ' » ‘ f A 


3-/^o 

z P//>° 


_ i 


>>//-; 2 ) = + 

'/• ' 4/^. 4-^V 


8 f/f. 


At low density, the function D is approximately: 


I 


T> 


dr 


~~ L k 


1 


/ 


Z5(3-l) (fyj 
Thus, all isentropes converge toward the value 


2fS-l) 


e-e 0 


l 


C 7 - 25 (5-l) 

Since, in the same limit 


as -fi- 
A 


-l 

2S(s- 1 ) 


( 49 ) 


(50) 


( 51 ) 


(52) 


it follows that also approaches zero along an isentrope. In fact, 

it is clear from Eq. (35) that approaches zero more rapidly 

than does 

The C , 5 form of the state equation derived in Eq. (36) is based 

on the linear approximation to the ^5 , ^ Hugoniot, and applies only 

when ^ . Thus, it must be viewed with suspicion when f*/f 0 

becomes even as low as unity, to say nothing of the limit 7 ^ — ► O . One 

/ * 

of the most serious deficiencies of Eq. (36) at low density is that negative 
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pressures are encountered along the isentropes. For example, when 
S = 2 , the ^ and Q coordinates of an isentrope are given by (defining 

p = - £ + + ° w <»> 


& 4 4 <5 


+ <"l - ¥¥'+¥ +o ^> (54) 


Figure 3 shows the low-density behavior, for S = 2 , of isentropes start- 
ing from density ratios at the shock of 1. 2, 1.4, and 1. 8. For the lower 
entropy states (lower values of ) the negative pressures are encoun- 

tered at higher density. 

Usually, the negative pressures themselves can be tolerated, espec- 
ially if they comprise only a small part of an otherwise high-pressure 
solution. Their effect is often more serious, however, in that they may 
lead to imaginary values of the sound speed. The square of the speed of 
sound, equal to the slope of the curves in Fig. 3, is negative over a rather 
wide range for the two lower-pressure isentropes shown. 

In the analysis below, Eq. (36) is used for all values of , in 

spite of its improper behavior at low density. The justification for use of 
such a model lies chiefly in the fact that the low-density states have no 
appreciable effect on quantities such as the shock trajectory. Furthermore, 
even if their effect were appreciable, it would be very difficult to select 
more accurate model. The states reached by a shocked material when 
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it is expanded to low pressure are understood far less fully than are the 
corresponding high-pressure states. If the material is completely vapor- 
ized, it may be possible to treat it as a perfect gas, but in other regimes, 
processes such as condensation render such an approximation invalid. 

It is reasonable to expect that the shock-propagation solutions are 
insensitive to the low-density portions of the equation of state only for the 
s emi-infinite -ta rget problem, because most of the mass is concentrated 
in the high-pressure region near the shock, where the equation of state is 
valid. Since only a very small fraction of the total mass is described by 
the low-density portion of the state equation, the over-all error in quanti- 
ties such as the shock trajectory is quite small. 

The situation is not nearly as favorable in the case of the thin-plate 
perforation problem. There, the momentum distribution of the material 
ejected from the rear surface of the plate is an item of central interest. 

In predicting such a distribution, the mass and velocity acquired upon ex- 
pansion to low pressure must be known accurately, and they may depend, 
in a sensitive manner, on the precise form used for the equation of state. 

It is noteworthy, in this connection, that most theoretical treatments of the 
thin-plate problem (see Ref. 19, for example) make use of the perfect-gas 
approximation. 

Negative pressures are present in the state equations used for some 

20 

of the computer solutions of hypervelocity impact. Riney, for example, 
extends the state equation into the tensile regime by use of the formula 

t- + (55) 
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where is the state equation used for ^ .> 1 . The equa- 

tion of state for iron used by Bjork^ also predicts negative pressures. 
Whether the same remark applies to his calculations in aluminum 1 is un- 
certain, since the state equation used in those calculations has apparently 
not been published. 

Walsh and his coworkers make use of a low-density equation pre- 
1 7 

sented by Tillotson in which the pressure is expressed as the sum of a 

perfect-gas term plus a correction which vanishes exponentially as the 

density goes to zero. The boundary along which the state equation is 

changed from the condensed-state form to the expanded- state form has 

4 

been slightly changed in recent work from that used in previous publica- 
3, 1 7 

tions. In both formulations, it is possible to encounter negative pres- 

2 1 

sures. Shanfield, Lee, and Bach, for example, have recently found 
negative pressures in aluminum when applying Tillotson's formulation. 

Figure 4 indicates the magnitude of the negative pressures that may 
be encountered. Shown here are the pressure-density states through 
which an iron particle passes after being compressed to twice normal 
density, according to the equations employed by Tillotson, Bjork, and 
Riney, and also according to the present formulation (Eq. (36)) using 
5 = 1. 59, = 1. 259 megabars. For convenience, the actual equa- 

tions used are collected in Appendix A. 

The fact that no difficulty caused by these negative pressures has 
ever been reported indicates that calculations in a semi-infinite target 
can safely be made. That such errors in the state equation can be toler- 
ated in the thin-plate problem is much less likely. 
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E. Sound Speed and Temperature 


As a final note on the thermodynamic properties of a C , S medium, 
it is of interest to calculate the sound speed and temperature field. The 
most convenient expression to use for the former is 



Use of Eq. (36) leads, after a little algebra, to the expression 


^ , (zs . 0 l/tll +. ^ (57) 

£ ///» 

The temperatures reached by a compressed solid have been calcu- 
lated in a number of different ways (see, for example, Refs. 22, 23, 24). 
Since a general form of the equation of state is being employed here, it 
is of interest to inquire whether the associated temperature field can be 
dete rmined with a comparable degree of generality. For this purpose, 
the method advanced by Walsh and Christian in Ref. 24 is especially con- 
venient. These authors determined the temperatures along the Hugoniot 
from the relation 


[ = e ' - e - + J ** 

° Muf (58) 

= ^ fc-O + f 

■r. n- 

where IT — Vq is the specific volume, and where the integration is 
carried out along the Hugoniot. On this path, the right-hand side of (58) 
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is known, and is, in fact, a function only of VJ - . Differentiating with 
respect to Vj 1 gives 


— f WoU} 

A i L ^ 






+ 

Z 


(59) 


An independent relation for the left-hand side of this equation can 
be found from 


Td^ = c v ctr + ( 2 ^ tVit (60) 

Integrating along the Hugoniot, and differentiating with respect to 
gives 



Equating (59) to (61) produces a linear differential equation for 7j" 




+ 



(62) 


where 





Assuming that 


C v and 



are constants, the solution is 


(63) 


\T 




t>V“ 


£55^1/“ 




(64) 


where the parameter b is related to the Gruneisen parameter at nor- 
mal density, i. e. , 
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Along the Hugoniot of a C , S medium 




(s- (s-. )/>/£,] 


( 65 ) 


( 66 ) 


A little algebra shows that 



z- S — ^ 5-H ) ^ 

[s- (s- 1 ^ 3 5 


• f(V) = - 5/> 0 


.x 

[s ~ ] 4 


(67) 


whe re 



Substituting these into Eq. (64) leads to 



whe re 


( 68 ) 



to 




<r(<r-»f dbr 


(69) 


It is clear from this analysis that the temperatures along the Hugoniot 

are controlled by the Gruneisen factor ^/p 9 and by characteristic 

temperature ^ • The values of these parameters given by Walsh 

24 

and Christian for Aluminum, Zinc, and Copper are listed in the table 
below. 
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Having the temperatures along the Hugoniot, the temperature at any 
other state can now be determined by finding how the temperature falls 
along an isentrope. This relationship has been given by a number of 
authors (see Ref. 16, for example), and is rederived in an alternate way 
in Appendix B. It has the form 

f - •* B Vj 

For the particular case of a C , -S medium, this takes the simple form 

f • (if" 

It is interesting that this expression has the same form as that of a per- 
fect gas of constant specific -heat ratio / 

f = (il <»> 

in spite of the presence of the term A 

F. Comparison with Enig's Formulation 
„ . 25 

Emg has recently derived an analytical formulation of the equa- 
tion of state of a C, 5 medium by pursuing the logical consequences of 
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the "mirror image" approximation. This approximation assumes that 
the pressures and particle speeds experienced by a particle in a planar 
isentropic expansion wave are mirror images of the Hugoniot curve: 



This approximation is consistent with another well-established approx- 
imation that the free -surface velocity acquired by isentropic expansion 
from a point on the Hugoniot is twice the particle 9peed on the Hugoniot. 
Enig noted that the "mirror-image" assumption implies an equation of 
state, and found a closed-form expression for this equation for the case 
of a C , S medium. In the present notation, Enig’s result is 


4s r 


e-eo 

c u 


whe re 





(73) 


* - ( l+ *&) <! ’+[ 4s ^' 1 )] 1741 

and where the functional dependence of ^ on is given parametri- 

cally by 


t 


(2i 


J. 


4e 
1+2 . 



(75) 
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This equation of state is shown in Fig. 5. It is noteworthy that only the 

parameters 5 t/^ 0 c‘ u ’ 3 (fe ~‘) ’ and appear. There 

is no other dependence on the parameter 

Enig goes on to calculate the temperature, using the assumption 
that the thermal coefficient of volume expansion at low pressure is con- 
stant. The resulting specific hear at constant pressure exhibits a rather 
singular behavior at low pressure, indicating that a satisfactory state 
equation cannot be derived under the simultaneous assumption of the mirror- 
image approximation, the linear U s , U t relation, and the constancy 
of the thermal coefficient cf volume expansion at low pressure. The 
formulation involving only g , ^ , and ^ is perfectly satisfactory, 

however, and can be used in its present form if one is willing to dis- 
pense with its temperature predictions at low pressure. 

The present form of the Mie-Gruneisen equation does not recover 
the mirror-image approximation. To illustrate this, calculations were 
made, for S = 2 , of the free-surface velocity acquired by a shocked 

particle. Along an isentrope, = - dClJjl , or 




c J A q. c 

° y c 

where subscript FS stands for free surface, and where the sound speed 
is given by Eq. (57). Because of the improper behavior of the present 


formulation at low pressure, the sound speed is imaginary over an appre- 


ciable range of pressure, as shown in the sketch below. 
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If the integration is started at the point A, the results are as follows 


% 

Vc 

a, /c 

Ufs /u, 

1. 8 

9 

4 

3. 04 

1. 4 

7 

2/3 

2. 10 


The principal contributions to the integral in Eq # (76) come from the 
low-density range. It is clear from these results that the present for- 
mulation would be quite inadequate for any study (such as the thin-plate 
problem) where low-density behavior is important. It also serves to 
indicate that caution must be used in interpreting results from computer 
solutions whose low-density state equations allow negative pressures. 
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HI- METHODS OF NONSIMILAR SOLUTION 

This section describes three approximate analytical methods that 
have been used for predicting some of the features of the flow patterns 
set in motion by hypervelocity impact. The objective of these solution 
methods is to produce simple formulae for the location and strength of the 
main shock wave, and for some of the details of the flow pattern behind 
the shock. None of the three methods uses a similarity assumption, and 
thus their success is measured by their ability to achieve maximum accu- 
racy with minimum complexity, by making the proper approximation. 

It should be emphasized that the only class of problem considered 
is that for which the projectile may be represented by a point disturbance. 
Thus, the analysis will not be valid for greatly elongated projectiles. 
Moreover, even fcr cubic or spherical projectiles, the analysis only 
applies when the scale of the disturbance is large compared with the scale 
of the projectile. 

A. Oshima's Method 

The literature of blast waves in gases contains a number of analytic 

methods for dealing with the nonsimilar effect. Notable among these are 

2 A 

perturbation methods, explored by Sakurai, among others. Such per- 
turbation methods, which describe the first-order departure from simi- 
larity, have been applied, for example by Kochina, 27 to shock propagation 
in condensed media. Another method, due to Oshima, 6 ' 28 offers the 
possibility of a wider range of validity than that of the perturbation tech- 
niques, for a comparable expenditure of effort. Furthermore, the method, 
developed originally to handle nonsimilarity arising from counterpressure 
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^ in the perfect-gas case, can easily be adapted to handle that result- 
ing from the equation of state in the present case. 

The first of the subsections that follow contains a review of 
Oshima's technique for the perfect-gas case, together with a few addi- 
tional features not explicitly mentioned in his original paper. The second 
section below then presents the application of Oshima’s technique to the 
case of shock propagation in a C , S medium. 

1 . Review of the Perfect-Gas Case 
a) General Analysis 

Consider the case of a symmetric flow, where the index is 0 

for the planar case, 1 for the cylindrical case, and 2 for the spherical 
case. Oshima begins by making the change of independent variables from 
f and £ to y) and . Thus the Euler equations become 


( 4 , - y ?) 




-h 




(77) 


* 2 - 




i 3 £ _ g £><£>_ 


(78) 
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These equations contain no approximations beyond those implied by use 
of an inviscid fluid. The terms on the right-hand sides come from the 
time derivatives in the original equations, and represent the change in 
the dependent variables with respect to the scale of the disturb- 

ance. For the case of a perfect gas, where 

? = m] 

the energy equation becomes 


& 


f + 2£l m i£ 1 


(81) 


The boundary conditions at the shock can be written conveniently in terms 
of the shock Mach number 


M = 


- Vc. 


(82) 


where the sound speed ahead of the shock, C Q , is 


= 


A 


(83) 


In terms of , the Rankine -Hugoniot conditions for a perfect gas are 





%s) 


2 *"-/ 

r+i y(r+i) m 7 - 

(V- »W**+z- 


( 84 ) 
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At this point, a similarity solution would require three steps. First, 

derivatives with respect to must be assumed to vanish. Second, 

the combination ^ mus ^ be assumed constant, which implies 

> fJ 

that l ‘ Finally, the shock Mach number must be taken to be 

large, so that M Z is small compared to unity. 

Oshima uses the notation 




__ &L _ /-A / 


/e 5 




(85) 


The value of fj is not assumed to be constant, however, as it would in 
a similarity solution, but is permitted to vary with (-L) . The 

next step is to approximate the derivatives with respect to /£ s . This 

is done by first replacing these derivatives by derivatives with respect 
to the shock Mach number, and by then approximating the result. If 
stands for either cf> , , or ^ , the final form is 


/? 22 = &32L fS. - . td k SZ = _ H I2C- 

3 24 24 U&, A/ ^ 24 A / 3rt 


( 86 ) 


using the definition given in Eq. (85). The problem is thus converted into 
one of approximating 3S /2M • (» is possible, of course, to proceed 

to this point directly from the basic partial differential equations, by 
using yj and as independent variables. ) 

The essential feature of Oshima's method is his approximation for 
, namely 


9* 


’Vs 


dH 


(87) 


Wj value, of 


V 


7|=, 
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The Rankine -Hugoniot conditions can be used to evaluate the right-hand 
side. Using (87), Oshima has, in general 



Note that the right-hand side depends only on 7^ , with as a param- 

eter. Thus the equations of motion have been reduced to a set of ordinary 
differential equations, which are found, after a little algebra, to be 


4 // ■+ 'P ( +- ->> y" ) = 


1-hJ -4 

V z+(Y-i)m z ' 


z /—A/ 


(89) 


(90) 


-2 


/-fs/ 

4 




/-^ jf zfjM 4X ) 

4 ^ (91> 


These equations differ from those of the similarity analysis by the 
presence of the parameter M • For co . they reduce to the 

self-similar relations which were first derived and solved numerically 
by G. I. Taylor. ^ 

Equations (89) to (91) can be solved for the derivatives in the form 


±3_ 

(4>-y\ V"- 


(92) 
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(93) 


4>'= 


- -^7-l. K - \p -*$ 




/ 

{ = 


(4>--nf f?s - +- ^ 

V%/p 


(94) 


whe re 


je = 


i-hJ 




P — l) 


P4> 


(95) 


1 V pi 1 - 1 * 


(96) 


& « 


/-V 

a/ 


z + 


2.^-1) 


4Y 


ZXt4 >r -(x- 1) z.4-(y'-0M 1 


£ 


(97) 


For a given value of ^ , these equations must be solved anew for 

each value of pi .In doing so, the value of a/ for each given value 
of y] must be chosen. Oshima points out in his original article that, 
for each ^ » there is only one value of a / which produces a solution 

free of singularities. When /s/ differs from this value, the solution 
ounters one of two singularities, because either or 

-Y[ vanishes. The existence of a single value of f\f leading to a 


enc 
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regular solution has been demonstrated both by Oshima's own calcula- 
tions for ~t) - 1 , ^ and also by some more accurate calculations of 

30 

Lewis for 7 ) =0, 1 , 2. Both these authors terminate the integration 
at a small, nonzero value of 7^ , which they refer to as the core 

radius. The value of which they find varies from the classical 

strong-blast limit of at M = <X> , to the value // = 1 at 

M = 1 , consistent with the fact that the shock expands like a strong 

blast wave at first, and ultimately decays to a sound wave, for which 
uo t . 

The distributions of pressure and density found by Oshima's method 
for 7 ^ = 2 bear a strong resemblance to the results of machine solu- 
tions of the full partial differential equations. Lewis' results, for ex- 
ample, can be compared with the correlated numerical solutions shown 
by Sedov (Ref. 15, pp. 246-248). The velocity distribution is not nearly 
as well recovered, however; the Oshima result has OL always positive, 
whereas the numerical result shows both positive and negative regions. 
The more gross properties of the solution, such as the predicted shock 
trajectory, agree well with machine solutions, and also with experiment. 

In order to determine the shock trajectory, the conservation of 
total energy is used 

A, N , 

E- - % j (e-e. 

(98) 



fry 4 ° 
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where 


and 



( 99 ) 


= 2 J 2 tt 4ir -for i> = ot,l. 


( 100 ) 


The constant £_ has the dimensions of energy, energy per unit length, 
or energy per unit area, according as - 7 * = 2 , 1 , or 0 . Equation ( 98 ) 
can be written in a form valid for any £ and 




4 


'*-0 


(V-#-0(V-0 
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( 101 ) 


where, as in all explosion problems, the scaling length depends on the 
energy release divided by the characteristic pressure in the medium 

\6>+r> 

Once the shock radius -shock speed relation is known, the shock trajectory 
is found from the identity 


F 0 = 






( 102 ) 



(103) 


b) Qualitative Nature of the Flow 

The role played by the singularities of Eqs. ( 92 ) to (94) and the 

general nature of the flow they define can be seen byrewriting the equa- 

1 5 

tions in the (X, Of ) coordinate system used by Sedov: 
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w04 



.= jC*£ 

v ^ 


In terms of these variables, the equations of motion can be written in the 
fo rm 


+ F ^_ -5. 


• . 0 :• 


<*¥- 






d.Jb*y^ 


[( x - a /) u - yy)fj - 6r-i)^g- + (r-Q«fr(<-^)FT 
(*v)V y 


( i 06 J 


( 1 0 7 * 


where the functions 


ff 


? and depend only on x and ^ 


F t = -(W/) x + 


4 (i — a /) 

2 + (Y_,)m x 


< 1 08 


V = - x (x -a/) + (l—d) 


mV/ 

mV 




(.09 






4Y 

z+(r-0^. 


2 (x- 



j 
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The use of the Sedov variables achieves the same result here as it does 
in the similarity limit, namely, the variable does not appear ex- 

plicitly. Thus, by dividing Eq. (107) by (106), a single differential equa- 
tion for V* as a function of ;< and ^ can be derived, and is useful 
in any detailed study of the singularities of this set of equations. The 
boundary values at the shock are 


5 Y+i 




r-i 


x+i (i r+C)Yrt\ 


i ■+■ 


(Y-f )H 


(in) 


The singularities of this system of equations occur at the points 

3 0 

y-zV =0 and ('y_/\/) 2 L tj, = 0. A plot of the solutions found by Lewis 
for a typical Mach number and for various values of has the appear- 

anc e: 



For values of a/ different from that which gives <£)( o ) = 0 , 

the solution encounters one of the singular lines, as shown in the sketch. 

The nature of the flow is the same as that of the similarity solu- 
tion 18 (fiSf=ac>) . The quantity t always remains negative, 

approaching — OO as 0. 
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c) A First Integral 

The similarity equations ( Yl-eO ) have a first integral, that has 

32 

been derived, among others, by Lees and Kubota. The quasi- 
similar equations have an analogous integral. It can be derived by writing 
the basic equations in the form 

V + ^ +- J < 112 > 




(113) 


('V - -Jj-r P') = ^ (114) 

where 

H, - * (ns) 

^ Z + tY-^M*- 




/-a/ mV/ 

a/ m 7 "- / 


(116) 




0 


4)f 


Z-h (Y’-i)M 1 


(117) 


If (1 1 2), divided 

Hjfrt-V) / . 

H.) 


by . is added to (114) divided by 

, the resulting equation can be integrated, to give 
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(118) 


^ -f* / 4-^1 l-ql 7 * = 


whe re 


<*= ('**-*•)/„ 


(119) 
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The same first integral has also been treated by Oshima, but his re- 

* 

suit appears to contain a misprint. 

34 

For the similarity case ( ) Von Neumann and, independently, 

1 5 

Sedov have derived an exact solution of the entire problem by express- 
ing, in terms of the similarity variables, the statement that the time rate 
of change of the energy enclosed between two surfaces must equal the 

rate at which work is being done on the surfaces, less the flux out. The 

35 

same exact solution has been later rediscovered by J. L. Taylor, 

▼ + , 36 „ , .37 

Latter, bakurai. 

For the quasi-similar case the same approach does not 

produce an exact solution, 
d) Mass Conservation 

Oshima's approximation does not conserve mass. The integrated 
value of the mass distribution from the origin to the shock should yield 
the same value as the product of the undisturbed density and the volume 
swept out by the shock, i. e. 


S 


I t 


f — T, 


A 




* The exponent on the quantity L (x-f") >c "l should be multiplied by CL 
in Eq. (24) of Ref. 33. V f > 
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J s 

{ ^ = 777 < 12 °> 

The alteration that is made to the continuity equation in the quasi-similar 
approximation violates this fact. The continuity relation, Eq. (112), can 
be written in the form 


4 (+*>)• '-T>^ + H4' 

Multiplying by -y ^ and integrating by parts gives 


(i 2i ; 


(V+l-H,) J ^ 4^ ( 4>~y) 

0 1 


( 122 ) 


If the limits are taken from zero to one, and the first of Eqs. (16) is 
used, the result is 


J — 


1 + 1 - H. 


(123) 


The right-hand side of this equation, evaluated from Lewis' calculations, 
is shown in Fig. 6 for 1 = 2. The maximum error, about 40%, occurs 
near - 2. Such an overestimate of the total mass present is quite 
significant, and is partly responsible for the distorted particle trajec- 
tories described below. 

e) Approximate Solution for yl y- 0 

The perfect-gas solutions reported by Oshima and by Lewis show 
a constant, nonzero pressure near the origin ( ^ = 0), while the density 
vanishes as a power of 
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} -f Oz. ~f 0 y '-p ~ ^ < 124 ) 

Substituting these into the continuity and energy equations ( ( 1 1 Z ) and 
(114)) gives, respectively, the relations 


(fl-l) b + f\ (n)+ /) — /■/, y — X t ) L — 


whose solution is 


4 = 


hf,+ H*/V 

7> + l 


7 >-M - H, - 


(125) 


(126) 


In order to completely specify the solution, it is also necessary to know 
and ^ • There appears to be no way of determining these analy- 

tically. They are related to each other, however, so an estimate of -f- 

o 

(perhaps from numerical solutions for a/ close to the proper value, and 
close to the origin) can be used to find an estimate of £3 . The 

relation, which comes from the first integral described above (118), is 


=[<Po))‘ ^ [{(>)}* 1 4>o)-ij 


(127) 


2. Asymmetric Disturbances in a Solid 
a) General Analysis 

The foregoing discussion deals exclusively with symmetric disturb- 
ances, such as a spherical blast. Because the impact problem is asym- 
metric with respect to the target surface, the £ -variations in the equations 
of motion must be reinstated before any meaningful results can be ex- 
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tracted. When such variations are permitted, however, it is necessary 
to solve a set of partial differential equations, even if the similarity 
assumption is made. At present, there are no analytic solutions account- 
ing for asymmetry, either with or without the similarity assumption. 

Limited analytical progress has been made in the similarity limit, 
however, by restricting attention to the axis of symmetry. 18 Approximate 
solutions along that axis have been made by approximating the off-axis 
pressure distribution. It was found that for a wide range of off-axis dis- 
tributions, the solution was virtually the same, and was quite close, in 
many respects, to the solution of the symmetric problem. That evidence, 
which was restricted to the perfect-gas (or extremely high pressure) 
limit, suggested that the impact process could be considered as one half 
of a spherically symmetric disturbance, whose scale is controlled by the 
energy of the projectile, its momentum playing a secondary role. Walsh 

and his coworkers have recently presented an exact numerical solution 

4 

of the perfect-gas problem. Several important features of their results 
are matched closely by the approximate solution. 

An approximate treatment of the asymmetric, perfect-gas problem 

38 

has also been given by Rayzer. His solution is quite close to the sym- 
metric, cons tant- ene rgy solution in many important respects. 

Finally, mention should be made of two other papers in which the 
asymmetric problem is discussed. In one of these, a brief reference 
is made (p. 63) to the use of a dipole representation. The other paper 41 

The case of impact into a half-space of gas represents the spherical 
generalization of the planar "sharp-blow" problem, which is discussed 
in detail by Mirels. 
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is chiefly a presentation of experimental data on the cylindrical problem 
of shock propagation due to a sudden energy release along a line (an ex- 
ploding wire) at the free surface of a half-space of water. The paper 

refers to an analytical treatment of the problem in a Russian disserta- 
42 

tion. 

The main conclusion drawn from the perfect-gas studies was that 
approximate solutions for quantities like the shock trajectory could be 
obtained by assuming the disturbance to be equivalent to one half of a 
spherically symmetric disturbance. It was emphasized, at the time that 
conclusion was presented, ^ that the evidence for it was limited to the 
perfect-gas, or similarity limit. In taking up the nonsimilar nature of 
the problem, then, one of the first questions to be answered is whether 
the same conclusion can be reached at finite values of the shock speed. 

The evidence presented below indicates an affirmative answer to this 
question. After the evidence has been reviewed, the discussion then 
centers on the symmetric solutions. 

The basic differential equations in the 7^ , & , J^ s coordinate 

system are given above in Eqs. (11) - (14). To these must be added the 
equation of state, which in these variables is 


r- 


-p 


'p-i 


2(S-I) I ^ 


(128) 


The boundary conditions at the shock for a C , S medium, written 
in terms of the shock Mach number M = ^ s / c , take the form: 


to, f ('(-4) 


(129) 
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W ^ — 

/ 4 - (s-l)l^f 


( 130 ) 




ZS 2 - 



(131) 


For a spherical shock, there is no tangential component of velocity imme- 
diately behind the shock. Thus 


= O 


(132) 


Using the same notation as for the perfect-gas case, a new symbol 
is introduced for the decay coefficient 


£ 


/-a/ 

a/ 


(133) 


and the Oshima approximation is again made, by use of the relation 




92 * 

9Z S 


/-'V 

rJ 


M 




&(') 


3 - 


(134) 


The derivatives o£ 4> , ^ . and f are readily iound from this 

formula. It cannot be used for the derivative of , however, since 

that function is identically zero along a hemispherical shock. This deriv 
ative is retained, for the moment, in its original form. 

There are two possible ways of calculating the derivative of ^ 

The first is to apply Eq. (134) to the boundary value of , giving 
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( 135 ) 




4 - 

M-l 


The equation of state, (128), can now be used to replace ^ in terms of 
-f and ^ . The seond way of calculating the derivative of ^ is to 

reverse these two steps, substituting (128) first, and then applying for- 
mula (134). The result is completely different 


^ - 'zli L ( I £ + 

s 3 %* 4 z(s-i) ( rf-i [s-(s- 1 ) 4 ) 


-f 


-b 


L ** [s- (s-l)^] 2 - 


/ + ( s-i)rt 


(136) 


These two approximations are referred to in what follows as the first and 
second forms of Oshima's method. The former has the advantage of 
greater simplicity, while the latter appears to be a more consistent for- 
mulation. Only the first form has been used in the calculations reported 
here. 

The same ambiguity is present in the perfect-gas case. ^ In the 
derivation outlined in Eqs. (79), (80), (81), and (91), the equivalent 
second form was derived. The first-form result is quite different, and 
in fact is considerably more complex. 

When the first form for a C , 5 medium is used, the basic equa- 
tions become 
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( 137 ) 
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^ M-l a/ / ±(s-}")rf 

For the case of a symmetric solution, where all £ -derivatives 
or for a "centerline" solution (i. e. , along the axis 9 = 0 in the 
trie case) these take the form 


(4>-n) P' + P ( V +■ V ~ + -£■ 


/ -<y 4p_ 

a / t + (S-l)M 


(138) 


(139) 


(140) 

are zero 
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(141) 
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( 142 ) 


^ 4^ d M-l 


. z ^^-±(^y. z *JL- -if- 


- a / -f /' P 


(143) 


(s-i)H 


where tr is defined for p = 2 as 

2uJ 


2r - Z 


9<9 


(^ o) 


(144) 


A different definition would be appropriate for a centerline solution in the 
cylindrical case ( 7 ) = 1 ), while for a planar case, *£T must be identically 


zero. 


If the function is now expressed in terms of -f and ^ , 

energy equation is 

1-rJ H _ t-J 


the 


p'-A z 4^' = 


A/ M'l 






(145) 


where the function A , defined as 


= — = +tx 





^ [s-(s-l) 4/\ *" 


(146) 


is a dimensionless sound speed (squared), i. e. , (see Eq. (57)) 


47 


AI-1821-A-2 



( 147 ) 



With these definitions, the final form of the qua si-* similar equations is 




(148) 


(4-n) V + -4 = ^ 

T 


(149) 







(150) 


where 



I- rJ / 

1+ (s-i)H 
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Solved for the derivatives, these equations become 


, , (4>-y)Kt -h Z-z -'I'fZL 

- 





( 154 ) 


(155) 


(4-v 


(156) 


These equations are identical in form to those which applyfor the perfect- 

1 8 

gas case, to which they reduce as - » CO . They have the same 
singularities as those for M — t 22 ; one is at the point where <$> =y\ 

and the other is at the ,! sonic line”, ^ where K-'ii is equal to the local 
sound speed. 

The Oshima approximation for a C , S medium fails to conserve 
the total mass, just as in the perfect-gas case. Following the same 
series of steps as is used in deriving Eq. (123), the result (which applies 
only for a symmetric case) is 


whe re 



/-A/ / 

^ i + (s-t)H 


(157) 


(158) 
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b) Results of Calculations 

Equations (154) to (156) were solved numerically on an IBM 704 
computer' for the spherical case ( V = 2), with 3 = 2, and using the 
first-form definition of /<^ . Spherically symmetric solutions ( •£* = 0) 

were done for M = 8, 4, 2. 5, 2, 1. 5, and 1.25, and "centerline” solu- 
tions ( “t" 0) for = 8, 4, 2, and 1. 5. The results are shown in Fig. 
7. The solution for M = CD would be the same as the perfect-gas 
similarity solution^ for = 3. 

The solutions with V ^ 0 require an additional equation specify- 
i n § 'C • In the absence of an exact solution, some approximation must 
be used. The approximation used here is to relate 'C to the off-axis 
pressure distribution in the same manner as used in the similarity solu- 
tion -- i. e. , the term /^ s ^ in Eq. (139) is assumed to be 

zero, the resulting differential equation is differentiated with respect to 
& , and then 0 is set equal to zero. The result is 


' Z 


t—h! 

rV 


V -h 


i(*-n 0 ^— 

^ J 4.^ zy[ 


i 

n.4" ze*’ 


(H °) = 


<d 


(159) 


The pressure variation is then approximated in the form 


££ 

9^ 


* - 


( 160 ) 


The value /<( = 1 0 , which was found satisfactory in the perfect-gas 
case, has been used in all the calculations reported here. With this speci- 
fication, the variation of 7T can be found as part of the solution. 


J# 

The author is very grateful to Mr. Harold M. Rosenbaum, Computer 
services Department, for his very competent handling of these calculations. 
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Near the shock, the only important terms in Eq. (159) are 


-k 






(’'*) f 


( 161 ) 


whose approximate solution 


tr 


i<£0) 0-n) 

[ i - 


= Kfo) ii-n ) 2 


(162) 


explains the parabolic variation of TT that is evident in Fig. 7. 

For each value of pS] , various values of a/ must be tried until 
a solution is found which is free of singularities. In general, only one 
value of hi leads to such a solution. When hj is less than this num- 
ber, the quantity ($-"*{) —ft 1 * vanishes at some nonzero value of 7^ , 
and when hi is greater than this value, the solution encounters the 
singularity at 4> = y • 

The results shown in Fig. 7 suggest that, just as in the similarity 
( [^1 = CD ) limit, the three-dimensional effects are not large, at least as 
far as quantities like the shock trajectory are concerned. The evidence 
for this conclusion comes from two sources. First, as will be seen 
below, the shock trajectory is mostly affected by the mas s -ave raged 
values of velocity, pressure, and density. Because such averaging intro- 
duces a factor ^ into the integrands, only the values very close to 

the shock are influential. Thus, the relatively large differences in the 
distributions of <p , , and -f between symmetric ( *£? = 0) and 

asymmetric ( 0) solutions are not likely to be reflected in comparably 

large differences in shock trajectory. 
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The second line of evidence concerning the insensitivity of shock 
trajectory to asymmetry lies in the values of a/ . Regarding hJ as 
a function of , the shock trajectory can be found by integrating Eq. 

(85) in the form 

d %r — F(m) d63) 

where (— = , and where the normalizing constant is 

defined below. The fact is that a/ has approximately the same value, at 
given f^l , for both the symmetric and asymmetric solutions. Thus the 
shock trajectory is not markedly affected by symmetry considerations. 

It is interesting to observe that the value of a/ for the asymmetric 
solutions is always less than its corresponding value for the symmetric 
solutions. At the perfect-gas limit, Walsh 4 has pointed out that the value 
~ • 375 a 8 rees with the result CL = -j~j - . 59 found from his 
calculations. It appears that the present results for exhibit the same 

ratio as that for = CO . 
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fdasyn 1 
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. 375 

. 667 

. 600 

. 600 
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. 444 
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1. 447 
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The fact that this ratio is approximately constant is presumably related 

to the constancy of (X revealed by Walsh's solutions. There is at 

❖ 

present no analytical proof of this assertion, however, 
c) Further Discussion of Symmetric Solutions 

Given that the symmetric solutions can be used for a crude approx- 
imation to the shock trajectory, it is of interest to examine in detail their 
distributions of pressure, density, and particle velocity. These results 
are shown in Fig. 7. As f^J gets closer to one, a narrow region appears, 
in which large gradients of velocity, pressure, and density occur. These 
large gradients are not as sociated with singularities of the differential 
equations (the solutions are free of such singularities), and the question 

39 

arises whether they represent the appearance of an entropy line or 
perhaps of a second shock. In connection with the latter possibility, it 
must be recognized that, in a gross sense, the terms which are intro- 

43 

duced by Oshima's method are not unlike the artificial viscosity terms 
which are often used for the express purpose of generating shock waves. 

In order to identify the source of these steep-gradient regions, a 
calculation was made of the trajectories followed by various particles. 
There are two ways of calculating such trajectories. One is to integrate 
the relation , while the second is to follow lines of constant 

entropy (since the entropy of a particle does not change, once the particle 

5,c It should also be observed that was found to increase with 

increasing Y > and hence with increasing 3 , in the perfect-gas 

( = CO ) results of Ref. 18. The same tendency presumably persists at 

finite values of M’ , indicating that materials with larger values of 5 
will be associated with larger values of a/ and OC 
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has crossed the shock). For an exact solution of the partial differential 
equations, of course, both methods would lead to the same trajectories. 
Within the framework of Oshima's approximation, however, they do not 
necessarily lead to the same result, and the closeness between the two 
may be used to judge the validity of the approximation. 

The first trajectory calculated is that of the shock. Conservation 
of energy gives 


0 / 

= T(m;s) 


(164) 


whe re 


and 




(165) 



-for -J = /, 2 


(166) 


^ is defined here as one-half the value used in Eq. (100), to facil- 
itate application to impact problems, where it will be assumed (as in 
Ref. 5) that all the available energy goes into one half of a symmetric 
disturbance. It is interesting to note how the quantity d’ 2 ' presents 
itself as the proper quantity to use as an ambient pressure. This pres- 
sure, and the total energy £. , can be scaled out of the problem by the 

definition 



( 167 ) 
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In these terms, the shock radius -shock speed relation is 


-g-= /U * V) (168) 

5 

Figure 8 shows the result for S = 2 , compared with the earlier 

Quasi-Steady solution. The shock trajectory itself, found by integrating 
the identity 


C~t 

Zo 



(169) 


is shown in Fig. 9. 

This figure also shows particle paths, found as lines of constant 
entropy. This is done by first finding from Eq. (46) the variation, with 
-yj , of the quantity ^ / c?~ * for a given value of , i. e. , at a given 

,A i 

instant. From Eq. (48), it is also possible to find the values of C 
that exist along the shock. A graphical solution of these two relations 
then makes it possible to identify, for each at a given , the 

value 0 which the shock Mach number had at the instant when that 

particle was traversed by the shock. The results are shown in Fig. 10. 
Note that horizontal lines (constant t^\ 0 ) trace the history of a given 

particle. Thus, the particle intercepted by the shock ( = 1 ) at the 

instant when = 2. 5 has dropped back to the point = . 75 at the 
later instant when f^j = 2. Using these results, one can find the coordin- 
ates f/{£ = n. ^7/?. , » along a particle path. The results, 

shown in Fig. 9, suggest that the steep-gradient regions are associated 
with an entropy line, and not a second shock. 
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While the trajectories of the high-entropy particles (those processed 
by the strong portions of the shock) are reasonable, those crossing the 
shock at later time make no sense. All particles should initially follow 
the shock -- i. e. , they should be deflected to the right. At late time 
^ f*"1 however, the particles appear to start in the opposite direc- 

tion. In addition to this obvious fault, these trajectories do not agree 
with what would be had by integrating the velocity field. In all the solu- 
tions of Fig. 7, the quantity <£> is always positive. However, the 
constant-entropy lines of Fig. 9 show regions of negative velocity. 

In spite of its internally inconsistent description of particle trajec- 
tories, the Oshima solution makes an accurate prediction of the shock 
trajectory. The path traversed by the shock in Fig. 9 is, if anything, 
closer to experiment than the Quasi-Steady solution presented previously. 5 

In general, Oshima's method appears to work best at early time, 
when most of the mass is concentrated near the shock. During the weaker 
stages of shock propagation, the errors introduced by its approximations 
become more important, and its predictions of distributions of physical 
quantities less reliable. Throughout the entire range of shock strength, 
however, the integrated quantities required for shock-trajectory deter- 
minations are adequately handled. 

B. Quasi-Steady Method 

In the foregoing discussion of Oshima's method, the point has been 
made repeatedly that accurate shock-trajectory predictions require the 
distributions of pressure, density, and velocity only very near to the 
shock. Oshima's method is specifically designed to yield the exact values 
at the shock for these variables, and excellent approximations to their 


AI-1821 -A-2 


56 



first derivatives with respect to . Thus, any other approximation 

which is accurate at the shock is likely to be adequate for shock-trajectory 
calculations, even though its analytical basis might not be as rigorous as 
that of the Oshima method. One such approximation called the "Quasi- 
Steady" method, was advanced in Ref. 5. In the sections that follow, this 
method is reviewed, extended to later, weaker stages of shock propaga- 
tion, and compared extensively with experiment and with computer solu- 
tions. 


The Quasi-Steady method is based on the observation that, for 
shock-trajectory calculations, the only quantity needed is the integral X 
in Eq. (165), which depends on the shock speed with material prop- 

erties (such as S ) appearing as parameters. The Quasi-Steady approx- 
imation equates this integral to the value 2T, that it would have for a 
perfect-gas similarity solution with / so chosen as to give the correct 
values of pressure, density, and particle velocity at the shock. In the 
perfect-gas similarity solution, the boundary values at the shock are 
given in terms of j by 


*• - £ 


= -±- 

f>0<£ 


Z 

r-H 


> 


(170) 


A _ Tf+t e-e„ _ 2 . 

r-i ; u? (r+ o" 


The derivatives at 7? = 1 are not quite exact, since they depend on 
hj , which in turn depends on the entire approximation scheme over the 
whole range of 'ij from zero to one. 


It should be noted that the Quasi-Steady method cannot be used to 
describe the nonsimilarity arising from counterpressure in the perfect- 
gas problem. The reason is that the Quasi-Steady method always has 
4>= -f at the shock, whereas they are in general different for a 
perfect gas (see Eq. (84)). 
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In applying the Quasi-Steady solution, the Hugoniot values (JL X , pi * 
t are regarded as known, and the corresponding value of 
is then found from any one of the relations 

4 - +/ , 

f= 2. £-1 = Z -I - -y ~ -i (171) 

For example, when the shock speed is such as to create a density ratio 


/ ' j P Q =1.5 , the integral X~ is assumed to be the same as it would 

1 _ I.T+t 

be in a similarity solution with =. j-p " - 5“ . At a later instant 

when J ^ = 1* ^ , Y is taken to be 6, etc. For a point-source solu- 

tion, the shock speed is initially infinite, and the density ratio f*\ / 
is equal to its limiting value. Thus ^ starts out with the value 

As the shock speed decreases, yA j ^ 
decreases, and thus y increases, approaching infinity as ^ ap- 
proaches one. The large values of / that are thus encountered 


in an 


attempt to account for nonsimilarity constitute, in a sense, the counter- 

/ 44 

part of the strong-blast theory in which $ is allowed to approach 1. 

It should be emphasized that the Quasi-Steady method can be quickly 
applied to any material whose Hugoniot is known. A graph, for example, 
° f /i vs. Us will give y vs. iyL ^ , from Eq. (171). The corresponding 
values of X f , found from Fig. 11, are then used in Eq. (164) to give 
the shock radius -shock speed relation 


'/a 




Z-rf'ki x; [r(&)) 


(172) 
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The shock trajectory is then found from the identity 


r 

i “ / 


^4 


( 173 ) 


o 

These relationships hold for any Hugoniot curve. 

For the special case of a C , S medium, the shock-trajectory 
relations can be worked out once and for all in general form. The basic 
equations are the same as Eqs. (167) - (169): 


3_ _ 

& 




M 


4 / ^T = / 


(174) 


whe re 


xfn; s) = T,[^(ms)J 

and where is defined by 


(175) 


/= £2£i !h±tl d76) 

M-l 

Curves of Y[^\ an( ^ vs * //P f ^ or S = 1. 2, 1. 5, and 2. 0, 

were presented in Ref, 5. Those results were based on values of 
extending up to 20, and provided results out to around /g = 7 , 

where 1/^ was between . 8 and . 9, and /g 0 was approximately 4. 
The calculations have now been extended to weaker shock strengths, by 
adding the cases = 40, 70, and 100, for which the distributions behind 
the shock are shown in Fig. 12. The new shock-trajectory results, 
shown as the solid lines in Figs. 13 and 14, extend to ^ » 

where li es between . 96 and . 98, and where C ^/y^ is about 15. In 

addition, curves have been prepared which show the pressure and density 
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at the shock as a function of the depth to which the shock has penetrated. 
These results, which can be found from Eqs. (174), (38), and (40), are 
shown in Figs. 15 and 16. For convenience in re-plotting, coordinates 
of all these curves are given in Tables I and II. 

It is also possible to find the distributions of various quantities 
behind the shock, using the distributions for the various values of / 
given in Figs. 7a - j of Ref. 18, and in Fig. 12 of the present report. 

For example, the pressure distribution can be given as a function of •''/£? 
and by using the relations 


X = mVA?) 


f __ 


(177) 


± m 

These calculations have been done for Y Y 20, and for 5 =1 2, 1.5, 
and 2. 0. The results are shown in Figs. 1 7a - c. Unfortunately, there 
are very few numerical solutions available at present with which to com- 
pare these results. 

Solution at Early Time 

At early time, a point-source solution must contain a finite amount 
of energy in an infinitesimal volume. Thus the initial shock speed is 
predicted to be infinite. In actuality, the disturbance in the near vicinity 
of the impact point is not well represented by a point source, but at first 
exhibits a structure that is governed by the details of the geometrical 
configuration in which the projectile and target come together. Shortly 
after contact is made, the shock entering the target tends to become 
planar and parallel to the target surface. Its plane portion is soon 
affected by rarefaction waves originating at the intersection of the pro- 
jectile with the free surface, but the interaction near the impact point is 


AI-1821-A-2 


60 



roughly one -dimensional until the time when the shock has propagated 
into the target a distance on the order of the projectile length. During 
this initial period an exact solution for the (constant) shock speed can be 
found, from the Hugoniot curves of the projectile and target (see, for 
example, Ref. 18, page 17). The Quasi-Steady theo ry is improved by 
using the constant- speed solution up to the point where it intersects the 
point-source curve, using that curve thereafter. Curves showing con- 
stant multiples of the stress -wave speed, designated » have been 

added to Fig. 13, which presents shock trajectories for the special case 
of a C , S medium. 

Solution at Intermediate Time 

Once the early-time phase has been completed, the shock speed 
ceases to be constant, and makes a transition to the point-source curve. 

If the impact speed is sufficiently high, so that the early-time shock speed 
is many times the stress -wave speed, then the transition to the point- 
source solution occurs in the strong-shock portion of the point-source 
solution, where the shock radius expands in proportion to the 2/5 power 
of time. As the shock continues to expand, it slows down further, and 
the power of time according to which it expands (essentially the slope of 
the log-log £ £ curve) becomes larger, approaching the constant stress - 

wave speed ( cs?t) at late time. For a less severe impact, the tran- 
sition to the point-source solution occurs at sufficiently late time that 

the 2/5-power-law variation is not detected. The speeds presently attain- 

. Vi- 
able in laboratory experiments are not sufficiently high to produce a c 

variation (even if the very early time data could be properly resolved). 
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Such variations do appear, however, in computer solutions. Bjork's data 

2 

for iron striking tuff at 30 km/sec, for example, show this behavior, 

as do the results of Walsh et al for an iron-iron impact at 40 km/sec. 

Thus, in general, the shock expands according to the first power 

of time, both early and late. During the intermediate period, a variable 

power of the time is observed. The power may initially be as low as 2/5 

for a very high-speed impact; it is higher for more moderate impacts, 

and in all cases tends toward 1 at late time. 

It should be pointed out that a different interpretation has been ad- 

45 

vanced by Heyda and Riney. These authors believe that the solution 

for the late stages of the hydrodynamic flow is of the form 

* 

with the transition to the form ^ = C occurring only when material- 
strength effects become important. The present results indicate that, 
even without any strength influence, the wave speed in the inviscid solu- 
tion asymptotically approaches the acoustic velocity of the fluid being 
considered. In support of the variation £ s co t. ^ , Heyda and Riney 
present a large number of computed results for Aluminum -Aluminum and 
Lead-Lead impacts. However, the data (which are listed in Appendix C) 
only extend to shock radii on the order of 5 . Reference to 

Fig. 13 will show that the slope of the log ^ , log •£ curve is indeed 
around 2/3 in this region. It appears that the calculations must be ex- 
tended to later times ( /? s > 10 ^ ) before truly asymptotic conditions are 
achieved. 

Solution at Late Time 

For the case of a point-source disturbance in a perfect gas, where 
nonsimilarity is caused by counterpressure, the solution at late time 
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finds the wave speed approaching the sound speed ( M 1)» while the 
excess pressure (shock pressure less the ambient pressure) decays in- 
versely as the first power of the shock radius. The present problem 
involves nonsimilarity due to the form of the state equation only; the 
counterpressure term is always neglected. Thus it is not surprising to 
find in the present results that, while 1 at late time, the law of 

pressure decay differs from the acoustic ( ^ 60 /fe ) relation. Instead, 
the pressure decays as approximately the -7/4 power of 

This result can be demonstrated by noting that the function JZ, (Y) 
for Y > 1° is closely approximated by the relation 


r, fir) = 



(178) 


where = 0. 366, H = 1.7. As Y\ is allowed to approach one, 'Y' 
becomes very large 


M 



r = 


Zs 

M-l 


Thus 


X, IrM] = (m-i)'’ 

(*V) 

The shock- radius, shock-speed relation then becomes 



or 


— /A 
v* ( z 0 


-% 


(179) 


(180) 


(181) 


(182) 
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This expression then gives the asymptotic law of pressure decay as 


fs tt = _ -z 

For /) = 1 . 7 , = 0. 366 , this becomes 

-1.74.T 


& 




( 183 ) 


4 * a /f 


/■ 


(184) 


It should be noted that if r(M) were proportional to ( M’ - 1)^ 
as M -** 1 » then the pressure would decay as the -3/2 power of 

It is quite possible that a theory more accurate than the Quasi-Steady 
approximation would indeed reveal such a variation. For example, it is 
reasonable to assume that, at late time, the quantity +• £ e£ Z is always 
equal to its value at the shock. Thus T~ would be given by 

and the -3/2 power-law decay of shock pressure results. The same re- 
sult would be found if _L, (y) were assumed inversely proportional to 
■ 


The present results are not capable of resolving exactly what the 
pressure-decay law is. However, it appears not to be the inverse first - 
power familiar in gaseous acoustics, and is probably closer to the 3/2 
power mentioned above. 

Comparison with Experiment and with Calculated Solutions 

Figure 13 shows a correlation of the experimental data available at 
present on the propagation of impact-generated shock waves. Included 
are some measurements taken from a set of high-speed photographs of a 
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Lucite target, made by Mr. Robert Piacesi of the Naval Ordnance Lab- 


oratory. * Also included are the numerical results of Bjork for iron 

2 1 

striking tuff at 30 km/sec and iron on iron at 5 km/sec, the results of 

3 45 

Walsh et al for iron-iron at 40 km/sec, and Riney’s data in Aluminum 

and Lead, referred to above. Particularly for the high-speed numerical 
results, the different regimes are clearly shown. Most of the experi- 
mental data lie in a range so close to the stress -wave speed that the shock 
trajectory is essentially a constant- speed line throughout. 

Data on shock pressure as a function of shock radius are shown in 
Fig. 15. One set consists of the results computed by Walsh and Tillot- 

son for iron striking iron at 40 km/ sec. Another set is a group of 

46 

experimental measurements made by Charest. In addition, Riney's 

45 

results in Aluminum and Lead are shown. The agreement is quite good 


in all cases, and indicates a decay law where 


csoJZ". for jfg 


1 


tending toward ^ ^ for on the order of ten * 

Charest 4 ^ notes that his data, as well as some calculations done 
for his case by Walsh and by Riney, are correlated by a -1. 6 power law, 
and is careful to stress that this law must not be extrapolated beyond the 
range where it has been observed. The correlation shown in Fig. 15 
indicates some guidelines for making such an extrapolation. If one moves 
in the direction of smaller 9 a more ra pi d decay is to be ex- 

pected. Very near the surface, of course, the pressure will level off at 
the value appropriate to the impact point. In the direction of larger * 


* These data have not been published previously. We are very grateful 
to Mr. Piacesi, and to Dr. V. C. D. Dawson of NOL for providing copies 
of the photographs, and for their permission to use the data in this 
report. 
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it is probable that the law is still valid, at least to the point 

where strength effects become important. Unfortunately, there is at 
present no reliable estimate of the point at which this occurs. 

The coordinate of all the points shown in Figs. 13 and 15 are col- 
lected for easy reference in Appendix C. 

Comparison with Shock-Speed Data in Wax Targets 
47 

Karpov has recently presented a refined set of wave-speed meas- 
urements in a wax target, and compares the data with the constant C , S 
form of the Quasi-Steady theory. Actually, the Hugoniot for the particular 
target material is not characterized by a single pair of parameters C and 
5 , but is better approximated by 

= 2. 91 + 1.48 0£ t (km / sec) for ^ ^ 20 kba r 

= 1. 85 + 7. 32 Li f (km/sec) for 20 kbar 

An improved Quasi-Steady prediction can be made by using the complete 
Hugoniot. The results are shown in Fig. 18. The agreement is improved, 
compared to that found using the constant values C = 2. 91 km/sec, and 
S = 1.48. The theory still overestimates the shock speed, however, by 
as much as 20%. 

It is possible that part of the discrepancy can be traced to strength 

# 50 

effects, which have been shown to influence plane -wave propagation. 

In the later portions of Karpov's data ( y&CvA) the density change across 

the shock is only a few percent. Thus a serious question exists as to 

whether or not an inviscid model has any relevance. 

Part is also due to the failure of energy scaling, and part may be due 
to a two-wave structure^, 49 Q f the main shock. 
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Summa ry 


On the whole, the Quasi-Steady theory makes a satisfactory pre- 
diction of shock-wave trajectories, and provides a quantitative estimate, 
in terms of lZ 0 , of the extent of the strong-shock, transitional, and 
stress -wave regimes. It can be quickly and easily applied to any mater- 
ial whose Hugoniot is known. Its chief drawbacks are that it is based 
on pure energy scaling, and that it gives no information about the three- 
dimensional character of the solution, although the distributions it 
predicts behind the shock appear reasonable near the shock and near the 
axis of symmetry. 

The deficiencies of the Quasi-Steady solution are essentially the 
same as those of Oshima’s method. Thus, in spite of its lack of rigor, 
the Quasi-Steady solution is probably preferable, because of its ease of 
application. 

C. Porzel-Zaker Method 

An approximate analytic solution for strong blast waves in gases 

7 

was developed some years ago by Porzel. His purpose was to account 

g 

for departures from ideal-gas behavior. In 1959, Zaker adapted the 
analysis so as to make it applicable to the nonsimilar problem, and pre- 
sented a solution for a point-source explosion in a solid. In view of the 
relative unimportance of asymmetry on quantities like the shock trajec- 
tory, Zaker's results also have application to the present problem. In 
the subsections that follow, the Porzel-Zaker method is reviewed, and 
modified slightly so as to incorporate the C , S state equation described 
above. 
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1 . Approximate Solution 


The basic idea behind the Porzel-Zaker method is the observation 
that the continuity and momentum equations involve only ^ , 

and IL . Thus, an approximation to any one of these can be used to 
find corresponding approximations to the other two. The basic quantity 
for which an approximation is chosen is the density, which is assumed 
to vary as a power of the distance from the shock: ' 



(185) 


where, as before, subscript 1 denotes conditions immediately behind 
the shock, and where the exponent , a function only of time, is 

found from the mass conservation law 


which gives 



(186) 


,I87 > 

This approximation reproduces other solutions to a remarkable degree. 
For example, the perfect-gas similarity solutions are extremely well 
recovered. In particular, the case Y = 7 , for which jA- =4/3 , 

A 

leads to t - 1 , which is the exact solution for that case. 

To find the velocity field, the conservation law is written for the 
mass lying between the origin and a given particle path, as shown in the 
following sketch. 

Solutions of this sort have also been considered briefly by Nowak. ^ 
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Here j r ^(b) is the Lagrangian coordinate of the particle which was ori- 
ginally situated a distance f 0 from the origin. Taking the time deriv- 
ative of this equation along a particle path gives an expression for 

/jyfc , which can then be equated to the value of given by 

the continuity equation. The result is 


= -/ jf - z *r 11881 


The solution of this equation for (X , using conditions at the shock to 
evaluate the constant, is then 


U 


= UL 



l- UJU -X 


(189) 


where 



^ ^ fit 


( 190 ) 
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The quantity Qi. is always negative or zero, since increases 

* • 

with increasing Q , while the factor is the decay 

coefficient (see Eq. (85)), which varies from -3/2 at infinite shock speed, 
to zero at — C. . For a C , -S medium, Eq. (190) takes the form 


oc (t) = 




(191! 


/e s (*‘i)[t + (z- omJ 

Zaker notes that oi is small in the strong-shock region /cK ~ — 5- 

and neglects it for that reason. The additional error introduced by this 
step during the weaker stages has not been studied. 

The next step is to find the pressure distribution, from the momen- 
tum equation 


/ Dt 

Using Eq. (189) for u , it is possible to show that 


22 / 

/e. 


Integration of Eq. (192) then gives the pressure 




(192) 


(193) 


[44 dr = f ■ f ^ ~ - jV ^ 


(194) 
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then leads, after 


Use of the above expressions for 
some manipulation, to 


^s/ 5 1 



f> and Vu -/jyt 





t 


4 Z 



«v 

/?. 


frr 

r (u v*s) 

U] I 

- 


z. / A r /fe » 

( £- +z 





By setting = 1 , this gives a formula for the pressure at the ori- 
gin in terms of conditions at the shock 


ft (t) - 

^s/*l 



(196) 



~TT 







o(^A 


2. Energy Field and Shock Trajectory 

All that is now required to complete the solution of the problem is 
a determination of the internal energy Q -£ 0 as a function of position 
and time. Once this has been done, the conservation of total energy can 
then be used to derive a differential equation for the shock trajectory. 

To find the internal energy, Zaker uses the Hugoniot curve to calculate, 
for any A> and £> , the ’’waste heat", or energy unavailable for 
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mechanical work. For the present formulation, it appears simpler to 
use the state equation derived above, since it gives & - e D explicitly in 
and ^ 

To illustrate how this can be done, and to show the type of differ- 
ential equation for the shock trajectory that results, consider the case 
where DC is neglected. The conservation of total energy, using the 
present C , 3 form of the state equation, can be written 


terms of 


t 


£ = t _±_ _ g«/,d 

Zir I?/ I (_ 2. -5-1 zs-2. 



(197) 


where 


fr 




jr. 


('zs-z) [s- 



£ = 




a, ^ 




>e. 


(198) 


Carrying out the integration leads to 




= f(£ s ; s) + <5 (4 ; s) e, dg. 


( 199 ) 


whe re 








Us~ 1) Z (f+S) 25-2 J [s- (S-i)^ j 


I 

1 






(200) 


- G,(/e s j s)(u,ik s -u.f) 




/ 
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The last term in Eq. (199) can be rewritten as 


• 2 . 


r? ip d-tf _ p C_ 

s dJb = cfrt *Lt - ** S T " "s" 


(201) 


Thus, Eq. (199) can be solved, so as to give an explicit expression for 
as a function of and ■ After integration, the 

5jC 

shock trajectory is then found by applying Eq. (173). 

## # 

If (PC is not neglected, the third-order derivative will 

appear -- for example, from the term . However, the differ- 

ential equation for total energy conservation never contains the independent 
variable t explicitly. Thus, the complete solution simply requires a 
number of quadratures. 

The integrations indicated above have not been carried out. The 
analysis given is intended only to show the contribution that an explicit 
form of the state equation can make to Zaker's solution. It would be very 
interesting to complete these integrations, and to compare the results 
more extensively with the Quasi-Steady theory. Presumably, the shock 
trajectory would not be markedly different, but the distributions of , 

yO , and 0L behind the shock would undoubtedly be better approxima- 
tions of the exact numerical solutions. 


* All of these steps could, of course, be put in dimensionless form, 
using only the variables * an( ^ /fP 

** One of the points that needs further study is the early-time solution. 
Zaker uses initial conditions related to an equivalent sphere^of water, 
and finds early radius -time variations of the form & % cs>t » where a/ 
is approximately 2/8 to 2/9, instead of the value 2/5 that would be ex- 
pected for a strong blast. 
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IV SHOCK PROPAGATION IN POROUS TARGETS 

The potential utility of porous substances for achieving improved 

resistance to meteoroid puncture has often been cited -- for example, 

in Refs. 18 and 52. Recent experimental findings in substances such as 
53, 54 55 

foamed plastics, ’ and fiber metals tend to substantiate the ex- 
pected improvement. One of the factors responsible for the improvement 
is that porous materials absorb more energy per unit mass than do solid 
materials. Thus they offer the possibility of more rapid attenuation of 

t -gene rate d shock waves. In this section, the Quasi-Steady method 
is used to derive a solution for the propagation of spherical waves in a 
porous medium. The purpose is partly to illustrate how easily this 
method can be applied to a somewhat unconventional situation, and partly 

to shed light on the degree of attenuation that may be achieved by poros- 

❖ 

ity. 

A. Hugoniot of a Porous Material 

The equation of state for the nonporous form of a given material 
can be used to predict the Hugoniot for any degree of porosity. This fact 
is pointed out in three recent Soviet publications, 57-59 an d details are 
worked out in those references for specific state equations. For present 
purposes, it is useful to find the Hugoniot for substances whose equation 
of state is the form introduced above (Eq. (36)). This equation applies 
for any degree of porosity; the fact that the constants C and 5 were 
determined from experiments on the nonporous form is irrelevant. The 

Zaker applied his theoretical approach some time ago,® and con- 
cluded that porosity was effective in achieving a more rapid attenuation. 
Some aspects of plane -wave propagation in porous solids have been 
studied recently by Payton. 56 
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constants, however, have no direct relationship to the shock-wave proper- 
ties of the porous substance. For example, c will not in general be 
the weak-wave speed. 


The initial density of the porous material is designated by ^/p 
and the porosity is defined as 


yy\ 


-f 

A 


> i 


( 202 ) 


The Rankine-Hugoniot equations for the porous solid are 


V e o- z(f> P /> H ) > fn (203) 

Eliminating the energy between the first of these and the state equation 
gives, for the Hugoniot, 




[s-fs-i s - (s-i hA//i] 


(204) 


With this relation between and , the other Rankine-Hugoniot 

equations can then be used to find U s and 0i i . Figure 19 shows 
typical results, for S =1.5 , and three values of the porosity param- 
eter, YA =1,2 and 4. To provide a common standard of reference, 
all quantities are made dimensionless with respect to the values y4> and 
C Several features are worthy of notice. The first is that a weak 
wave ( 6 ► 0) compresses the material back to its normal density, con- 
sistent with the fact that a porous solid has very little compressive strength, 
and is returned to its solid form by even a slight disturbance. The most 
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striking feature, however, is the behavior of strong shocks. A limiting 
density ratio is achieved, equal to 



S 

S-l 


(205) 


If the initial porosity is greater than , even an infinitely 

strong shock will not compress the material back to its normal density, 
although an extremely weak wave will do so. The reason is that the tem- 
peratures achieved in a shocked porous medium increase with increasing 
porosity. At high enough porosity, the decrease in density associated 
with these high temperatures is sufficient to keep ^ less than /Q> 

C 7 C Q 

This phenomenon has been observed in the Soviet experiments, ’ and 
was also reported recently^ in this country. 

A detailed analysis of the weak-wave limit shows that the shock 
speed and particle speed are proportional to each other 




vh-i 

YV\ 




(206) 


and both are approaching zero. 

B. Solution for Shock Trajectory 

Having determined the Hugoniot for a porous C , S medium, it is 
now possible to find a solution for impact-generated shock propagation by 
using the Quasi-Steady solution. Conservation of energy in a spherical 
blast requires 

E=J (e + ^i/jz-rr^dr = — K 2 " ZT (rf J wi,s) (207) 
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whe re 


/W = 'Vc 


T = 


-J 


y 4 ; 3 = e 4" ; 4= A^ ( 

In terms of the usual scaling length, ^ 
comes 


2tt^C 


0 


■4 


( 208 ) 


, this be- 


^ _ f s) 

^ / ~>n 


-'/, 


(209) 


The integral X is now approximated by equating it to the value it would 
have for a perfect gas whose value of Y rnatches conditions at the 
shock 


/ = 


+\ 


A 


(210) 


As the shock speed varies from infinity to zero, increases from 

(yv\+\) / 

25-1 to ^ (yA * The resu ^ n S variation of shock speed with 

shock radius is shown in Fig. 20 for 5=1.5 , ~Y*\ =1,2 and 4. It 

should be noted that the wave speed at a given depth is reduced consider- 
ably by the porosity. 

Another integration gives the shock trajectory, from 


J M 


( 211 ) 
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which is shown in Fig. 21. The reference quantities C and depend 

only on the target material and impact energy; they are independent of 
the porosity. Thus the curves shown provide a direct measure of the in- 


creased attenuation caused by the porosity. 

One interesting feature of the curves for W\ 1 is that the solu- 

. l/g- __ 

tion has the form ^ c^o £ both at early and at late times. Thus it 
never departs by very much from the classical Taylor solution for a 
strong shock. Mathematically speaking, such behavior is due to the 
finite range of variation of . Whenever ^ is constant, the solutioi 

is of the form 


X(yO 


C£>ns7 


( 212 ) 


At early time, cc 2S — 1 , while at late time, * 


(wA -V L 




Physically speaking, both of these limits must be viewed with some reser- 


vation until experimental evidence becomes available. At early time, 
the pressures generated in a porous target are less than in a nonporous 
one; thus the point-source assumption of instant destruction of the pro- 
jectile is not as easy to justify. At the late-time limit, inviscid theory 
does not apply, because of the low pressures involved. 

In any event, the evidence suggests that a porous substance will 
decelerate an impact shock wave more rapidly than the nonporous form of 
the same material. Whether this can be a useful design concept in pro- 
tecting against meteoroid penetration depends upon a number of other 


systems considerations, of course, but it would appear that further study 
of the concept is warranted. In particular, studies of wave attenuation 
by internally instrumented targets would be very informative. 
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V CONCLUDING REMARKS 


Status of Analytical Theory 

The main purpose of the research described above has been to 
develop analytical methods for predicting the history of an impact- 
generated shock wave as it travels through a target. The problem is 
made complex by two factors: first, it is multidimensional, and second, 

the state equation cf solids renders similarity methods inapplicable. 

The present results indicate that the multidimensional nature of 
the problem can safely be ignored, if attention is restricted to quantities 
such as the shock trajectory. This conclusion is based upon approx- 
imate solutions along the axis of symmetry; and their comparison with 
solutions which assume the impact-generated flow field to be equivalent 

to one half of a symmetric disturbance. Such a conclusion had been 

5 

reached previouslyby Rae and Kirchner in the limit of extremely high 
impact speed; where similarity applies; the present results extend this 
conclusion to the range where ncnsimilar effects are felt. 

To account for the nonsimilarity of the problem, three approximate 
approaches have been studied and extended. One of these, called the 
Quasi-Steady method, has been developed more completely than the other 
two. The result is a simplified prediction, based on energy scaling, from 
which the history of the shock location, shock speed, shock pressure, 
etc. , can be found for any material whose Hugoniot is known. This theory 
has proven successful in correlating all cf the measured and calculated 
shock-trajectory data available at present. The only remaining discre- 
pancies occur either at very early time, due to the plane-wave nature of 
the problem in the near vicinity of the impact point, or at very late time, 
due possibly to the onset of mate rial -strength effects. 
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One of the most important results of the Quasi-Steady theory has 
been the identification of the scaling parameter R. 0 , proportional to 

the cube root of the projectile kinetic energy divided by the character- 
istic pressure f> 0 C ** of the target, and the demonstration, by com- 
parison with experiment and with computed results, of the ranges of 
this variable in which various types of shock motion can be observed. 
Thus, the perfect-gas, or similarity solution, in which ^ , 

can be detected only for times less than 10 ^ ^°/c > while the acoustic, 

or stress -wave, limit is encountered at times greater than 10 
The intermediate, or transitional, regime is the one in which most 
present-day results lie. In all three regimes, the shock trajectory begins 
with a constant- spee d, plane-wave phase, and approaches the point- 
source solution only after the projectile has been destroyed. 

Suggestions for Further Research 

Future improvements in the state of analytical theory could most 
profitably be sought in two areas. The first is concerned with accounting 
f c r the two-dimensional nature of the flow field. When approximate 
descriptions of this feature of the flow are obtained, they will permit the 
restriction to pure energy scaling to be removed. 

The second area in need of improvement is the late-time portion of 
the problem. Even in the inviscid approximation, the law of decay of the 
shock pressure should be firmly established. Looking beyond that stage, 
the problem of material-strength effects is amenable to analytical study. 

A great deal of computational effort is presently being directed at this 
aspect of hypervelocity impact, and such analytical developments would 
be helpful in interpreting the numerical results. 
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It should be reiterated that the solutions described above are 


capable of supplying* with very little effort, a considerable amount of 
information -- for example, concerning the shock trajectory. Their 
advantage, in this regard, lies in their ability to handle unconventional 
problems very quickly. A significant example is the treatment, given 
above, of shock propagation in porous solids. Another such example is 
the application of the Quasi-Steady theory to the problem of shock propa- 
gation due to irradiation of a solid by an intense laser beam. 

Implications for Damage Effects 

No report on the mechanics of hypervelocity impact should close 
without indicating how the studies in question contribute to our under- 
standing of impact damage. While shock-propagation phenomena are of 
great interest in themselves, nonetheless the subject of most intense 
practical concern is the configuration in which the target material finally 
comes to rest. 

At least three types of damage may be considered, depending on the 
thickness of the target. For very thin plates, whose role is to spread 
out the impulse of the projectile over a wider area, the essence of the 
problem is to determine the momentum distribution of the material coming 
from the back of the plate. Although this case has not been specifically 
investigated here, nevertheless the present studies have called attention 
to the significant effect of the low-density portions of the state equation. 

For somewhat thicker plates, the most serious mode of deformation 
is that of spall fracture. Here, the present analyses supply important 
information, since the shape and amplitude of the incident compression 
wave must be known before the occurrence of such fractures can be pre- 
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dieted. The companion publication to this one discusses in greater 
detail the application of the present results to the spallation problem. 

Finally, for targets which are effectively semi-infinite in extent, it 
is desired to know the crater depth. Unfortunately, no theory based on 
the invisc id-fluid model can offer any direct evidence on this question. 
Physically speaking, the crater boundary takes shape along that line 
where the material strength predominates to such an extent over normal 
and shear forces that further removal of material ceases. This occurs 
long after the point where the strengthless model is no longer valid. 

To remedy the situation, some auxiliary criterion for defining the 
crater size must be chosen, and its choice is crucial. The current dis- 
agreement over the scaling laws for crater size^ 3 can be traced almost 
entirely to differences in crater-formation criteria, 5 since all of the 
published shock trajectories comply with essentially the same scaling law. 

The two most prominent crater-formation criteria in use are those pro- 
1 3 

posed by Bjork and by Walsh, and the present studies bear, indirectly, 
on both of these. 

Bjork identifies the crater boundary as the locus of points where 
the pressure is extremely low, and the particle velocities randomly 
oriented. The existence, noted above, of negative pressures in the state 
equation used by Bjork for iron raises a question about the validity of 
the low-pressure portions of his results. Use of the inviscid model at 

It should also be noted that the Hugoniot curve for iron given by Bjork 
(Figs. 3 and 4 of Ref. 1) is in error. It should be compared with the 
experimental results listed by Rice, McQueen, and Walsh (Ref. 16) and 
by Altshuler et al (Ref. 64) or with the Hugoniot determined from the state 
equation used by Bjork (Fig. 1 of Ref. 18). 
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low pressure is always a questionable procedure, and the uncertainty 
is compounded when the state equation itself is suspect. What is needed 
is a quantitative statement of the pressure levels involved, and of the 

extent to which they are affected by the state equation. 

3 

Walsh bases his criterion on the late-stage equivalence which he 
observes. His argument is that* if two different impacts give rise to 
the same flow patterns during the late stages of the inviscid phase, then 
the subsequent deformation during the strength-affected phases must 
also be the same. Thus, crater size is scaled by the same parameters 
as those which govern the late-stage equivalence. 

By comparison, Walsh’s criterion appears to the present author to 
be preferable. However, certain features of the present studies suggest 
that further examination of this criterion may be profitable. In particu- 
lar, much of the evidence on which late-stage equivalence is based 
comes from integrated quantities, such as the total fluid momentum. It 
has been observed above that such integrals (in the case of the energy 
distribution) are rather insensitive to conditions at points whose dis- 
tances from the origin are less than about 1 / 3 of the shock radius. But 
these are the particles that will eventually form the crater boundary. It 
would be a more sensitive test of late-stage equivalence if some param- 
eter could be chosen that does not heavily weight the influence of particles 
near the shock* but rather emphasizes the role of the relatively less 
dense material. It has been pointed out, in this connection (Ref. 5, page 
43) that the time scales for shock propagation in different targets under 
identical impact conditions appear to be different from the time scales 
that govern the ejection of material from the crater. 


83 


AI-1821-A-2 



In summary, it can be said that the approximate analytic solutions 
discussed above properly account for nonsimilarity, accurately predict 
the shock trajectory, and focus attention on the need for better informa- 
tion on the state equation at low pressure. Future improvements should 
seek a better description of the two-dimensional nature of the solution, 
and should begin to account for the influence of material-strength effects. 
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APPENDIX A 


STATE EQUATIONS 

In Ref. 2, Bjork presents the state equation which he uses for iron, 
crediting it to Dr. F. H. Harlow of the Los Alamos Scientific Labor- 
atory. It has the form 


— | l,l3o XIO 4 ^-!) 4- 31, IS X 10* (4^ - | )|^->| 


+ iTLa (tf'-i) + 443.4 (A-i) 


whe re 


4 -///>. 


0. 3?<54 +0. 5-3 c ?4 




+ /> 0 


and where the megagram, millisecond, meter system 


3 3 

of units is used. Thus the density ^ is 7. 86 gm/cm = 7. 86 Mgm/m . 

2 2 

The internal energy bas the units (m/msec) , or (km/sec) , 

2 -2 / 

and the pressure is in Mgm/(msec) -m or 10 Mbars. For 1 » 

20 * 

Riney uses the same relation, and also credits it to Harlow. For 

1 , Bjork apparently continues to use Eq. (A-l), while Riney's 

formulation -- Eq. (55) -- takes the form 


f 


' ,c Riney's version contains several factors involving parameters £ 
and . Dr. Riney has informed the author in a private communica- 

tion (September 15, 1964) that these parameters have been set equal to 
zero in all calculations to date. It should also be noted that the absolute 
value sign in the first line of Eq. (A-l) is absent from Riney's formula. 
For 4 ^ 1 , its absence makes no difference. 


= 1 1,180 *io 4 +-/> 0 (V-eo )[<?r?. / -f-/s ~ 4>8 

4 /C(e-e 0 ) u [ o.*9fa +o.S4o6>(4'-i)] — - - 

1 ) 9oo-h/> 0 (e-e^) 


(A-2) 
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The state equation used in Walsh's calculations has the form, for 


condensed states 




+■ 


e-ea 


-*-i4 +- 2. ('P-i) 


-t-i 


while for expanded states it has the form 


(A- 3) 


f = ^(e-e 0 ) 



In these relations, the pressure is in megabars, the internal energy in 

3 2 7 

Mbar-cm /gm, and the constants, for iron, are 


a = o. 5 

b = 1.5 

,4 = 1. 279 Mbar 
B =1.05 Mbar 

= 0. 095 Mbar-cm^/gm 

(X = fi = 5 
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APPENDIX B 


ALTERNATE DERIVATION OF TEMPERATURE 
ON AN ISENTROPE 


This Appendix contains an alternate derivation of Eq. (70). The 
method is due to Sedov (Ref. 15) who points out (p. 233) that because 


cLa 


de + Q/f) 

~r 


(B - 1 ) 


is a perfect differential, it follows that a partial differential equation 
for “p can be derived from the exactness condition. In general, if 


is exact, then 


d LL — Mdx + 


N,C f 


ihL - = 0 

If the state equation is given in the form 
then 


(B-2) 


(B-3) 


(B-4) 


dA- 



(B-5) 


and the exactness condition, after multiplication by t 1 " > has the 

form 



95 


AI-1821-A-2 



-^■his partial differential equation may be written in the form of three 
ordinary differential equations 


dt = £ (b-7 

r t-rfr), s(m, 

by using the method of Lagrange (see, for example. Ref. 65, p. 250). 
The method states that, in general, the problem of solving 

p ^P)(5r) + Qb-'i i<*){0-) - ^3-*) <B ' 8: 

6 o 

is equivalent to the solution of the system of ordinary differential equa- 
tions 



If this set has the two solutions 


(B-9) 



and 


then the general solution of (B-8) is 



(B - 1 0) 


U — -f(\r j 

where -f is an arbitrary function. 

For the Mie-Gruneisen equation 


e - 


A- 



(B - 1 1 ) 


(B-12) 
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one of Eqs. (B-7) is 


cLt __ rfr)4/> 


(B-l 3) 


T 


whose solution is 




(B - 14) 


The second of Eqs. (B-7) 


tztlMk 


^ = i ' '~C'f (B-15) 

^ Z 5 ” Of), 

is a relation for the v coordinate of an isentrope, as can be seen by- 
comparison with Eq. (56). Thus the solution of (B-15) is found to be 




4 t-L 


(B-16) 


which is the same equation for an isentrope already derived in Eq. (44). 
The general solution for the temperature field is therefore 


f J 


T = 

where A* is the entropy. In particular, along an isentrope, 


(B - 1 7) 


= e+b [- f A 

-r- (. f / 




(B - 1 8) 


which is the result presented above (Eq. (70)). 
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APPENDIX C 


COLLECTION OF SHOCK-WAVE TRAJECTORY DATA 

The data shown in Figs. 13 and 15 come from a variety of sources. 
This Appendix lists the measurements used, and the quantities employed 
in reducing the data. 

Q 

Eichelberger and Gehring present data for a Lucite target struck by 
a steel projectile. The data listed in Table III below were taken from 
Fig. 5 of their paper. 

The second set of Lucite experiments was reported by Halperson 
and Hall. 10 They state that their results have the form ^ = 8 
where ^ is in mm/psec, and K s is in millimeters. Integration 
leads to 

= (5. 62 X 10 5 -£ ) 4 / 5 (C - 1 ) 

where is in cm, and t in seconds. This formula has been used 

to calculate the ^ , t entries in Table IV. 

The third set of measurements in Lucite was taken by Piacesi (see 
p. 65). The projectile was a 1 /4-inch-diameter aluminum sphere. The 
measurements are listed in Table V. 

The first set of numerical results was reported by Bjork. 1 The 
quantities listed below as are actually the distance along the axis of 

symmetry from the impact point to the shock. Because the shock itself 
is spread out over several mesh sizes by the artificial viscosity tech- 
nique, the precise location of the shock is uncertain by several percent. 

The numbers quoted below in Table VI lie about halfway between the points 
of minimum and maximum pressure. 
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The same comments apply to the data from Bjork's calculation of 
2 

iron striking tuff. In reducing these data, the value C = 2.21 km/sec 
was used. This was found from the equation of state given by Brode and 
Bjork. ^ At low pressure, it has the form 


t * 




(C-2) 


where, if units are used, the constant has the units 5. 30 x 10‘ 


dynes j gm 


c m 


/EEL\ 

\ erg J 


1/2 


Equating (C-2) and Eq. (8) gives 




Z-^o L^g 


(C-3) 


Thus, using Eq. (7) 




(C-4) 


or 


as 


A 

r 


A - A (t£) 

Z/3> V° ' 




£* - 




(C-5) 


Using ^ =5. 30x10 SIB — and A* = 

sec -cm^ 

km/sec. The data is given in Table VII. 


(C-6) 


3 

1. 7 gm/cm gives C = 2. 21 
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The third computed trajectory is that for iron striking iron at 
40 km/sec, given by Walsh and Tillotson. The /<? s , £ points, given 
in J.able VIII, were found from Fig. 12 of Ref. 3, which presents the 
pressure ^ | behind the shock as a function of • These pressures 

were converted to shock speeds by use of the Hugoniot curve (Fig. 8 of 
Ref. 17). The resulting graph of */f^ vs was then integrated 

{by Simpson's rule) to give £? s vs £ . The shock speed was taxen to 

be constant, at its irnpact-pcint value of 3. 05 x 10^ cm/sec for 
0 < ^ ”S 1 - 36 cm. The resulting trajectory is listed in Table VIZI. 

Another set of calculated trajectory data has been given by Lake 
and Todd (Fig. 7 of Ref. 67). These calculations assume a spherical 
wave, and neglect the Q -variation. The results, given in Table IX, 
are for an iron projectile of mass 10 gm, striking aluminum at 4, 75 
km/sec (15, 600 ft/sec) but reveal an extremely rapid rate of shock propa- 
gation. For example, at the latest time shown (6. 9 x 10~ 10 sec) the 

average shock Mach number M = ^ S /ct is 82 - 7 * These data have 

d 

net been included on Fig. 13. 

46 

Charest has recently presented a set of measurements of the shock 
pressures developed in an aluminum target when struck at 7. 32 km/sec 
by 0.476 cm diameter aluminum spheres. His data have been reduced, 
and appear in Table X. 

45 

r inally, Heyda and R:ney have presented computed results for a 
series of aluminum-aluminum and lead-lead interactions. Measurements 


The & -variation is included in the formulation of Ref. 68, but the 
continuity equation used there is in error. The equation in question -- 

Eq (1) of Ref. 68 -- is presented for the case of cylindrical, not spherical 
polar coordinates (cf Eq. (1) above). 
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taken from their figures, and the associated values of 






are given in Tables XI through XVI. 
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Table I 

QUASI-STEADY PREDICTION OF SHOCK RADIUS, SHOCK-SPEED RELATION 


i/m r r 8 /r 0 M/v 2 


S = 1.2 .1 1.667 

.2 2 

.3 2.113 

.4 3 

.5 3.8 

.6 5 

.7 7 

.8 I I 

.886 20 

.9415 40 

.9662 70 

.9762 100 

S = 1.5 .1 2.33 

.2 2.75 

.3 3.29 

.4 4 

.5 5 

.6 6.5 

.7 9 

.8 14 

.858 20 

.926 40 

.958 70 

.970 100 

S = 2 .1 3.44 

.2 4 

.3 4.714 

. 4 5.667 

. 5 7 

.6 9 

.7 12.333 

.8 19 

.8 1 20 

.9023 40 

.9436 70 

.9604 100 


.347 

75 

4 

.636 

16.7 

3 

.959 

6.49 

2.4 

1 . 34 

3 . 12 

2 

1.813 

1.67 

1.714 

2.43 

0.925 

1.5 

3 . 28 

.511 

1 . 333 

4.667 

. 260 

1.2 

6.96 

. 121 

1 . 105 

10.75 

.0550 

1.052 

15 . 10 

.0301 

1.029 

18.72 

.0208 

1.020 


.449 

60 

2.50 

.795 

18.25 

2 . 14 

1 . 174 

5.20 

1.87 

1.619 

2 . 50 

1.667 

2 . 16 

t . 333 

1.50 

2.85 

.734 

1.36 

3.80 

.409 

1 . 25 

5 . 34 

.208 

1 . 154 

6.80 

. 130 

1 . 105 

10.62 

.0569 

1.052 

15.0 

.0307 

t .029 

18.65 

.0211 

1.020 

.581 

45.0 

1.82 

1.017 

10 

1.667 

1.475 

3.89 

1.54 

2 

1.88 

1.428 

2.62 

1.0 

1.333 

3.44 

.556 

1.25 

4.55 

.306 

1 . 178 

6.35 

. 156 

1 . Hi 

6.56 

. 145 

1.105 

10.45 

. 0600 

1.052 

14.85 

.0317 

1.029 

18.50 

.0214 

1.020 
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Table II 

QUASI -STEADY SOLUTION FOR SHOCK TRAJECTORY 


M R o 

ct / R 0 ) S =1.2 

ct / R o)s = 1.5 

ct / R o)$ = 2 

. 5 

.0356 

.0206 

.0198 

1 

. 153 

. 1128 

.0887 

1.5 

.342 

.2706 

.2144 

2 

.586 

.4826 

.3908 

2.5 

.873 

.7386 

.61 16 

3 

1. 194 

1.0319 

.8703 

3.5 

1 - 543 

1 . 3546 

1.160 

4 

1.914 

1.7016 

1.476 

4.5 

2.302 

2.0686 

1.814 

5 

2.704 

2. 4526 

2. 171 

5.5 

3. 119 

2.8496 

2.543 

6 

3.542 

3. 2576 

2.927 

6.5 

3.975 

3.6756 

3.324 

7 

4.415 

4. 1026 

3.732 

8 

5.312 

4.974 

4.571 

9 

6.225 

5.867 

5.437 

10 

7. 151 

6.778 

6.324 

1 1 

8.086 

7.702 

7.226 

12 

9.029 

8.635 

8. 140 

13 

9.981 

9.577 

9.064 

14 

10.939 

10.526 

9.997 

15 

1 1.901 

1 1 . 480 

10.938 

16 

12.868 

12.440 

1 1.882 

17 

13.839 

13.403 

12.833 

18 

14.812 

14.369 

13.789 

19 

15.788 

15.339 

14.734 
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Table J 1 1 

El CHELBERGER AND GEHRING (REF. 9) IRON - LUCITE 


M = 0.18 gm 

V = 4. 6 km/ sec 

c = 2.59 km/ sec 


R 0 = 0.338 cm 

R 0 /c = 1 . 304/z sec 

p 0 c i2 = 7.92 x I0 10 

dynes/cm 

tjjX sec 

R s , cm 

ct/R 0 

R c /R n 

S 0 

.6711 

.334 

.516 

.988 

1.06 

.540 

.613 

1.598 

1.55 

.77 

1. 19 

2. 27 

2, 10 

.92 

1.62 

2.72 

2.40 

1.07 

1.84 

3. 16 

3.20 

1.40 

2.45 

4. 15 

4. 30 

1.59 

3.30 

4.70 

5.6 

1.93 

4.30 

5.70 

6.9 

2. 23 

5. 30 

6.60 

8. 2 

2.56 

6. 30 

7.60 


HALPERSON 

Table IX 
AND HALL (REF. 10) 

ALUMINUM - LUCITE 


M = 1 .0 gm 
R 0 = .610 cm 

V = 5.31 km/sec 
R 0 /c = 2.35ytisec 

c =2.59 km/sec 
P Q c 2 = 7.92 x 10 10 
dynes/cm^ 

t, jj. sec 

R s , cm 

ct/R 0 

M R o 

1.780 

1.780 x 10 1 

1 

6.3 

0.755 

7.55 

1.640 

10.3 
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Table I 

PI ACES I (PRIVATE COMMUNICATION) ALUMINUM -LUCITE 

M = 0.361 gm V = 5.7 km/sec c = 2.59 km/sec 
R 0 = 0.1*90 cm R 0 /c = 1.90/tsec /) Q c 2 = 7.92 x I0 10 dyne$/cm 2 


t ,ji sec 

Rg t cm 

ct/R 0 

R A 

1.076 

0.75 

0.565 

1.53 

2.552 

1.50 

1.345 

3.06 

4.028 

2.05 

2, II 

4. 19 

5.504 

2.55 

2.90 

5.20 

6.980 

3. 17 

3.66 

6.48 

8.456 

3.69 

4.45 

7.51 

9.932 

4.20 

5.22 

8.59 

11.41 

4.71 

6.00 

9.61 

12.88 

5.08 

6.80 

10.3 


Table H 


' 

BJORK (REF. 1) 

IRON - IRON 


M = 6 160 gm 

V = 5.5 km/sec c 

= 4.0 km/sec 

R 0 = 4.90 cm 

R 0 /c = I2.25^i.sec 

2 = 1.259 Mbar 

t ,ji sec 

R s , cm 

ct/R 0 

R s/ R o 

3.5 

3. 1 

.286 

.632 

8.7 

7.9 

.710 

1.61 

81.7 

36 

6.68 

7.35 


Table VI 1 



BJORK (REF. 2) 

IRON - TUFF 


M = 1.065 x 10 

10 gm V = 30 km/sec 

c = 2.21 km/sec 

R 0 = ^6.5 meters 

R 0 /c = 2. 10 x 10 

sec /V 2 = 8 

.30 x I0 10 dynes/cm 

t, $ec 

R s , meters 

ct/R 0 

M R o 

. 17 x I0' 3 

5.5 

8. 10 x I0 -3 

. 1 18 

.36 x I0" 3 

11 

1.72 x I0‘ 2 

.236 

3.44 x I0" 3 

65 

1.64 x I0‘ 1 

1.40 

6.36 x I0" 3 

89 

3.02 x 10" 1 

1.92 

9.25 x I0" 3 

105 

4.40 x 10" 1 

2.25 

24.8 x I0* 3 

165 

1.18 

3.55 

61 x I0' 3 

260 

2.90 

5.60 


AI-1821-A-2 


106 



Table HU 

WALSH AND TILLOTSON (REF. 3) IRON - IRON 

M = 38.6 gm V = 40 km/sec c = 4.0 km/sec 

R 0 = 3.4 cm R 0 /c = 8.5/tsec /^o° 2 = 1-259 Mbar 


t.yu, sec 


R s , cm 


P i . Mbar 

ct/R 0 

R s/ R o 


.450 


1.37 



.053 

.403 


.678 


2 


- 

.0798 

.589 

- 

.894 


2.5 


- 

. 105 

.735 

- 

1.150 


3 


- 

. 135 

.881 

- 

1.449 


3.5 


- 

. 170 

! . 03 

- 

1.787 


4 



.210 

1. 18 

- 

2. 160 


4.5 


- 

.255 

1.32 

- 

2.566 


5 



.303 

1.47 

- 

3.012 


5.5 


- 

.355 

1.62 

- 

3.497 


6 


- 

.410 

1.76 

- 

4.019 


6.5 


- 

.474 

t .9 ! 

- 

4.586 1 


7 



.540 

2.05 

- 

5. 186 


7.5 


- 

.610 

2.20 

- 

5.812 


8 


- 

.685 

2.35 

- 

- 


1.34 


47.5 

- 

.394 

37.8 

- 


2.7 


17.5 

- 

.795 

13.9 

- 


3.8 


8.4 

- 

1.12 

6.67 

- 


6.2 


2.9 

- 

1.82 

2.30 

- 


6.5 


2.6 

- 

1.91 

2.06 

- 


7.25 


2.0 

- 

2. 13 

1.59 

- 


8.2 


1.5 

- 

2.41 

1.19 





Table IT 





LAKE AND TODD (REF. 

67) IRON - ALUMINUM 


M 

= 

1 Or 3 g m 

V = 

= 4.75 km/sec c = 

5.85 km/sec 


R o 


2.695 x 

I0‘ 2 

cm 

R 0 /c = 4.61 

x IO' 8 sec 


t,/*sec 

R s- 

cm 

ct/R 0 

R s/ R o 


CM 

in 

X 

.O'" 

9. 12 

x i<r 2 

1. 13 x IO" 3 

3.39 


1.6 

X 

10-'° 

1.52 

x 10' 1 

3.47 x IO' 3 

5.65 


3.9 

X 

IO-'0 

2.41 

X 10' 1 

8.46 x IO" 3 

8.95 


6.9 

X 

I0‘ 10 

3.35 

X 10* 1 

1.5 x IO" 2 

) 2 • 4 



107 


AI-1821 -A-2 




Table X 

CKAREST (REF. 46) ALUMINUM - ALUMINUM 


M = 

0. 1525 gm 

V = 7.32 km/sec 

C = 

5.85 km/sec 

*o 

= 0.187 cm 

Ro^c = 

0. 32(tyL sec 

/>o c2 = 

0.924 Mbar 

R g , cm 

p | , Mbar 

R s/ R o 



0.51 


.352 

2.73 


.381 

0.76 


.206 

4.05 


.224 

1.02 


.150 

5.45 


. 162 

1.02 


. 130 

5.45 


. 141 

1 . 53 


.0786 

8.20 


.0850 

1.91 


.0496 

10.2 


.0538 

2.54 


.0315 

13.6 


.0341 

3.05 


.0239 

16.3 


.0259 

3.55 


.0198 

19.0 


.0215 

4.80 


.0120 

25.6 


.0130 

6.40 


.0076 

34. 1 


.00822 

7.61 


.0053 

40.8 


.00573 



Table H 




HEYDA AND 

RINEY (REF. 

45) ALUMINUM 

- ALUMINUM 

M = 0 

. 1525 gm 

V = 7.6 km/sec 

c = 5.85 km/ sec 

II 

o 

0. 197 cm 

R 0 /c = 0.337^/ sec / 7 0 

c 2 = 0.924 Mbar 

t'^lsec 

R g j cro 

P|, Mbar 

c t/R 0 

R ./ R O 

P|//V 2 

.340 

.315 

1.155 

1.01 

1.60 

1.25 

.380 

.345 

1. 1 10 

1. 13 

1.75 

1.20 

.420 

.380 

.865 

t .25 

1.93 

.935 

.510 

.410 

.695 

1.52 

2.09 

.750 

.750 

.615 

.540 

2.23 

3. 13 

.584 

.915 

.750 

.240 

2.71 

3.81 

.260 

.985 

.790 

.225 

2.92 

4.01 

.244 

1.050 

.850 

.226 

3. i 1 

4.32 

.245 

1. 180 

.920 

.210 

3.50 

4.68 

.227 

1.320 

.990 

.21 1 

3.91 

5.02 

.229 

o 

in 

s* 

i. no 

.110 

4.30 

5.63 

. 119 
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Table TTT 

HEY DA AND RINEY (REF. 45) ALUMINUM - ALUMINUM 

M = .1525 gm V = 20 km/sec c = 5.85 km/sec 

R 0 = .375 cm R 0 /c = .642/* sec J* 0 c 2 = .924 Mbar 


t sec 

R 8 , cm 

p lt Mbar 

ct/R 0 

R ./ R 0 

Pl./’o 

.200 

.320 


<311 

.854 


.210 

.350 

4.88 

.328 

.935 

5.29 

.232 

.380 

4.30 

.362 

1.015 

4.65 

.260 

.415 

3.95 

.405 

1. 1 1 

4.28 

.285 

.435 

3.60 

.445 

1.16 

3.90 

.375 

.555 

1.52 

.585 

1.48 

1.64 

.410 

.630 

1.53 

.640 

1.68 

1.66 

.505 

.695 

1.37 

.786 

1.86 

1.48 

.555 

.770 

- 

.865 

2.05 

- 

.632 

.825 

0.94 

.986 

2.20 

1.02 

.750 

.955 

0.73 

1.170 

2.55 

.790 

.860 

1.075 

0.78 

1.34 

2.87 

.843 

1 .045 

1. 1 10 

0.43 

1.64 

2.96 

.465 

1.555 

1.370 

0.52 

2.42 

3.65 

.561 

1.823 

1.770 

0.23 

2.85 

4.73 

.250 

2.076 

2.025 

0.20 

3.22 

5.40 

.216 
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Table HU 


HEY DA AND RINEY (REF. 45) ALUMINUM - ALUMINUM 
M = .1525 gm V = 40 km/sec c 3 5.85 km/sec 

R 0 = .595 cm R 0 /c = 1.020 sec f> 0 c 2 = .924 Mbar 


t.^usec 

R s , cm 

c t/R 0 

R ,/ R 0 

• 122 

.315 

.120 

.530 

. 133 

.340 

. 130 

.570 

. 140 

.370 

. 137 

.622 

.156 

.410 

. 155 

.690 

. 165 

. 44 

. 162 

.740 

. 18 i 

.465 

. 178 

.783 

. 19 1 

.500 

. 187 

.840 

.252 

.62 

.247 

1.04 

.280 

.69 

.275 

i . 16 

.310 

.74 

.305 

1 . 24 

.366 

.82 

.360 

1.38 

.396 

.87 

.390 

1.46 

.43 

.92 

.421 

1.5 4 

.46 

1.0 

.470 

1.68 

.52 

1.05 

■ 510 

1.76 

.58 

1 . 15 

.570 

1.94 

.74 

1.25 

.725 

2 . 10 

.68 

1.35 

.863 

2.27 

t .03 

1.5 

1.01 

2.52 

1.60 

2 . 10 

1.57 

3.54 

1.74 

2.20 

1.70 

3.70 

1.90 

2.40 

1.86 

4.04 
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Table UZ 

HEYDA AND RINEY (REF. 45) ALUMINUM - ALUMINUM 


M = . 1 525 gm 
R 0 = .780 cm 


V = 60 km/sec c = 5.85 km/sec 

R 0 /c = 1 . 332/4 sec /> 0 c 2 = .924 Mbar 


R s , cm 

t.^ec 

ct/R c 

V R o 

.31 1 

.087 

.0653 

.400 

. 340 

.090 

.0675 

•435 

.377 

. 102 

.0766 

.485 

.141 

. i to 

.0826 

.526 

. 1*45 

. 12 1 

.0910 

.570 

.48 

. 131 

.0962 

.615 

.51 

. 141 

. 106 

.652 

.63 

. 180 

. 135 

.809 

.70 

.202 

.152 

.897 

.76 

.230 

. 173 

.972 

.83 

.260 

. 195 

1.06 

.90 

.282 

.212 

t . 16 

.95 

.331 

.249 

1.22 

1.0 

.370 

.278 

1.28 

1.15 

.44 

.330 

1.48 

1.26 

.56 

.420 

1.62 

1.40 

.63 

.474 

1.80 

1.50 

.70 

.526 

1.92 

1.67 

.84 

.630 

2. 14 

1.80 

.96 

.720 

2.30 

1.95 

1 . 1 1 

.832 

2.50 

2.08 

1.25 

.940 

2.66 

2 . 17 

1.40 

1.05 

2.79 

2.34 

1.52 

1.14 

3.00 

2.50 

1.74 

1.31 

3.20 

2.60 

1.90 

1 . 43 

3.34 
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Table JL 

HEYDA AND RINEY (REF. 45) LEAD-LEAD 
M = .641 gm V = 7.6 km/sec c = 2.07 km/sec 

R ° = -394 cm Rq/c = 1. 904 }l sec /° Q c 2 = .434 Mbar 


fJL sec 

R s , cm 

P|, Mbar 

Ct/R 0 

M R o 

P l ,f 0 c 

.40 

.285 

4.21 

.210 

.722 

8.70 

.45 

.310 

3.33 

.236 

.788 

6.88 

.50 

.335 

2.87 

.263 

.850 

5.92 

.55 

.375 

2.69 

.290 

.950 

5.56 

.60 

.405 

2.53 

.315 

1.03 

5.21 

.70 

.424 

1.78 

.367 

1.08 

3.67 

.80 

.485 

1.98 

.420 

1.23 

4.09 

.87 

.560 

1.23 

.456 

1.42 

2. 54 

.95 

.620 

.95 

.500 

t . 58 

1.96 

1.30 

.685 

.81 

.682 

1.74 

1.67 

1.62 

.885 

.61 

.850 

2.25 

1.26 

1.80 

.950 

.57 

.946 

2.41 

1. 18 

2.00 

i.OI 

.43 

1.05 

2.56 

.89 

2. 17 

1.07 

- 

1. 14 

2.72 

- 

2.31 

1 . 1 1 

.33 

1.22 

2.81 

.67 

2. 70 

1.24 

.37 

1.42 

3. 15 

.765 

4. 70 

1.76 

. 125 

2. 46 

4.46 

.258 

5.90 

2.03 

. 135 

3. 10 

5. 15 

.278 

6.29 

2.28 

.06 

3.30 

5.80 

. 124 
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Table XV I 

HEYDA AND RINEY (REF. 45) LEAD - LEAD 


M = 0.641 

gm 


V = 20 km/sec 


c = 2.07 

km/sec 

R 0 = .751 

cm 


R 0 /c = . 363yl£sec 

II 

CM 

U 

« 

.484 Mbar 

t ,/X sec 


R s , cm 

Pi, Mbar 

ct/R 0 

R ./ R 0 

Pi //> 0 c 2 

. 16 


.315 

21.2 

.0441 

.420 

43.9 

.24 


.350 

- 

.0660 

.465 

- 

.28 


.385 

16.3 

.0770 

.512 

33.7 

.30 


.410 

14.8 

.0829 

.545 

30.5 

.31 


.435 

13.85 

.0855 

.580 

28.6 

.36 


.495 

10.3 

.0992 

.660 

21.3 

.40 


. 560 

6. I 

. 1 10 

.745 

12.6 

.46 


.625 

5.3 

.127 

.830 

1 1.0 

.53 


.690 

5.45 

. 146 

.920 

11.2 

.60 


.755 

4. 1 

. 166 

1.00 

8.48 

.75 


.820 

3.6 

.207 

1.09 

7.42 

.90 


.950 

2.9 

.248 

1.26 

6.00 

1.05 


! . 08 

2.4 

.290 

1.44 

4.95 

1.27 


1. 12 

1.4 

•350 

1.49 

2.90 

t . 62 


1.25 

1 .7 

.446 

1.66 

3.51 

2.0 


1.50 

0.65 

.551 

2. 00 

1.34 

2.71 


1.77 

0.6 

.748 

2.34 

1.24 

3. 10 


1.895 

- 

.856 

2.52 

- 

3.37 


2.025 

0.45 

.930 

2.70 

.93 

4.20 


2.260 

0.40 

(. 16 

3.04 

.83 
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Figure I GENERAL SHOCK-WAVE RELATIONS FOR THE Ml E-GRUNEI SEN 
EQUATION OF STATE 
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Figure 5 ENIG'S EQUATION OF STATE 
















Figure 7d QUASI-SIMILAR BLAST WAVE SOLUTION 
>*= 2 S = 2 M = 2 
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Figure 8 SHOCK SPEED-SHOCK RADIUS RELATION FOR QUASI -STEADY SOLUTION 
AND QUASI-SIMILAR SOLUTION ; *> = 2, S = 2 
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Figure 13 SHOCK-TRAJECTORY DATA COMPARED WITH QUASI-STEADY SOLUTION 
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PRESSURE DISTRIBUTION ACCORDING 
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Figure 18 BLAST-WAVE THEORY COMPARED WITH NEW DATA IN WAX TARGETS 


INITIAL DENSITY OF NONPOROUS SOLID 







